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We have used two robotic telescopes to obtain time-series high-resolution optical echelle spectroscopy and VI and/or by 
photometry for a sample of 60 active stars, mostly binaries. Orbital solutions are presented for 26 double-lined systems 
and for 19 single-lined systems, seven of them for the first time but all of them with unprecedented phase coverage 
and accuracy. Eighteen systems turned out to be single stars. The total of 6,609 -R=55,000 echelle spectra are also used 
to systematically determine effective temperatures, gravities, metallicities, rotational velocities, lithium abundances and 
absolute Ha-core fluxes as a function of time. The photometry is used to infer unspotted brightness, V — I and/or b — y 
colors, spot-induced brightness amplitudes and precise rotation periods. An extra 22 radial-velocity standard stars were 
monitored throughout the science observations and yield a new barycentric zero point for our STELLA/SES robotic 
system. Our data are complemented by literature data and are used to determine rotation-temperature-activity relations 
for active binary components. We also relate lithium abundance to rotation and surface temperature. We find that 74% of 
all known rapidly-rotating active binary stars are synchronized and in circular orbits but 26% (61 systems) are rotating 
asynchronously of which half have P^ot > Porb and e > 0. Because rotational synchronization is predicted to occur 
before orbital circularization active binaries should undergo an extra spin-down besides tidal dissipation. We suspect this 
to be due to a magnetically channeled wind with its subsequent braking torque. We find a steep increase of rotation 
period with decreasing effective temperature for active stars, Prot oc , for both single and binaries, main sequence 
and evolved. For inactive, single giants with Prot > 100 d, the relation is much weaker, Prot oc T^ff''^^- Our data 
also indicate a period-activity relation for Ha of the form Rn^ oc Prot for binaries and Rua oc Prot*^* for singles. 
Its power-law difference is possibly significant. Lithium abundances in our (field-star) sample generally increase with 
effective temperature and are paralleled with an increase of the dispersion. The dispersion for binaries can be 1-2 orders 
of magnitude larger than for singles, peaking at an absolute spread of 3 orders of magnitude near Tcn ~5000 K. On 
average, binaries of comparable effective temperature appear to exhibit 0.25 dex less surface lithium than singles, as 
expected if the depletion mechanism is rotation dependent. We also find a trend of increased Li abundance with rotational 
period of form log n(Li) oc — 0.6 log Prot but again with a dispersion of as large as 3-4 orders of magnitude. 
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1 Introduction 



star arrives at the bottom of the red giant branch (Iben 1967 



In a previous paper (Strassmeier et al. 2000 paper I) we 
reported on radial and rotational velocities, chromospheric 
emission-line fluxes, lithium abundances, and rotation peri- 
ods of a total sample of 1,058 G5-K2 dwarfs, subgiants, and 
giants based on 1,429 moderate-resolution KPNO coude 
spectra and 8,038 Stromgren y photometric data points. The 
aim of this survey was to detect new candidates for Doppler 
imaging but, besides the discovery of 170 new variable stars 
and 36 new spectroscopic binaries, the more intriguing re- 
sult was that 74% of the G-K stars with Ca II H&K emission 
also showed significant lithium on their surface. However, 
G-K giants should have very few lithium on their surface 
because of convective mixing. Theoretical models predict 
that surface lithium has to be diluted by many factors once a 



* Based on data obtained with the STELLA robotic telescopes in Tener- 
ife, an AIP facility jointly operated with lAC, and the Automatic Photo- 
electric Telescopes in Arizona, jointly operated with Fairbom Observatory. 
** Corresponding author, e-mail: kstrassmeier@aip.de 



Charbonnel & Balachandran 20001. Out of the 21 Doppler 
imaging candidates found, just four stars were single stars, 
three of them evolved, the rest were spectroscopic binaries, 
but all four single stars had very strong lithium. 

Despite that it is generally acknowledged that higher 
than normal lithium abundance is common among magnet- 
ically active stars, no unique correlation with rotation rate 



was found after Skumanich's ( 1972 1 original discovery. Re- 
cently, White et al. ( |2007 1 revisited this issue in their sam- 
ple of solar-type dwarfs but no such correlation was found. 
Almost all surveys just revealed trends, if at all, and even 



these appear to be of different quality (e.g. Lebre et al. 2006 
Bohm-Vitense 2004; do Nascimento et al. 2000, 20(53; De 
Medeiros et al. ,2000, De Laverny et al. 2003 ; Randich et 
al. 1993[ l. The comprehensive survey of nearby giants by 



Luck & Heiter ([2007} did not even show a trend. How- 
ever, the line broadening in their stellar sample was just 
3-7 kms^^ and likely too narrow a range to see a trend. 



Bohm-Vitense ( |2004 1 suggested that the steep decrease of 
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V sin i in early G giants as well as in Hyades F dwarfs at 
effective temperatures cooler than «6,450 K, i.e. at their 
lithium dip at about the same temperature, are the result of 
deep mixing and related to the merging of the hydrogen and 
the helium convection zones. More recently, Takeda et al. 



(20101 announced evidence for a (positive) correlation of 
Li abundance with rotational velocity in a sample of solar- 
analog stars. 

The spectroscopic survey of 390 solar-like dwarf stars 
by White et al. ( 2007| l included 28 of the stars in our sur- 
vey. Relating Call H&K radiative losses to stellar rota- 
tion. White et al. ( |2007| l found a saturation of chromo- 
spheric emission for rotational velocities above approx- 
imately 30 kms~^. In an earlier paper, Strassmeier et 
al. ( 1994| l verified that evolved stars obey qualitatively the 
same scaling of Ca ll-K-line flux with stellar rotational ve- 
locity or period as do main-sequence stars (see, e.g., Mama- 



jek & Hillenbrand 2008 , Pace & Pasquini,2004 or Pizzolato 



et al. 2003 for a summary). No qualitative difference was 
found between single evolved stars and their equally rapidly 
rotating counterparts in a spectroscopic binary. However, 
large scatter indicated that rotation might not be the only 
relevant parameter. Based on a sample of 22 intermediate- 
mass G and K giants in close binaries, Gondoin ( |2007 1 not 
only verified the rotational dependency of (coronal) X-ray 
surface flux but also found a dependency on surface gravity. 
Such a dependence could stem from the effect of gravity on 
coronal electron density and on the overall sizes of coronal 
loops. 

Massarotti et al. (2008]! reported rotational and radial 
velocities for 761 giants within 100 pc of the Sun. They 
found that all binaries in their sample with periods less than 
20 days have circular orbits while about half the orbits with 
periods between 20-100 days still showed significant eccen- 
tricity. They also found evidence that the rotational velocity 
of horizontal branch stars is larger than that of first-ascend 
giants by a few kms^^. Earlier, De Medeiros et al. (2002i 



presented a study of 134 late-type giants in spectroscopic bi- 
naries and found a considerable number of G-K giant stars 
with moderate to moderately-high rotation rates. These ro- 
tators have orbital periods shorter than 250 days and circular 
or nearly circular orbits and appear to be synchronized with 
the orbit. 

The present paper follows up on the newly identified 
spectroscopic binaries with active components from our pa- 
per I. Its direct aim is to determine their orbits on the basis of 
high-precision radial velocities and to separate their compo- 
nent's rotation and activity tracers along with other absolute 
astrophysical parameters. Only with precise stellar parame- 
ters can we directly compare binary components with sin- 
gle stars and then be aware of the spectrum contamination 
from unknown secondaries or even tertiary stars. We recall 
that an unknown continuum contribution from a secondary 
star impacts on the determination of the effective temper- 
ature, gravity etc. and could together drastically alter the 
derived lithium abundances and thereby mask any relation 



if present. In Sect. |2]we restate our sample selection criteria 
and give a summary of the target stars. In Sect. |3] we de- 
scribe the new observations and in Sect. |4] we derive basic 
quantities from the spectra and the light curves. These in- 
clude radial velocities, orbital parameters, rotational veloc- 
ities and photometric periods, stellar atmospheric parame- 
ters like temperature, gravity and metallicity, lithium abun- 
dances, and absolute Ha-core fluxes. Sect. |5] lists notes to 
individual stars. Sect. [6] presents the analysis in terms of 
rotation, temperature, activity, and lithium-abundance rela- 
tions. Finally, Sect. |7] summarizes our findings and conclu- 
sions. 



2 Sample selection 

Our sample selection is based on the 1,058 stars from the 
KPNO Doppler-imaging candidate survey in paper I. It it- 
self was drawn from a total of 6,440 stars from the Hippar- 
cos catalog (ESA |1997} van Leeuwen |2007| l for the bright- 
ness range 7'."0-9™5 and declination -30° through H-70°, 
B — V colors between 0™67 and 1"^0 for stars with paral- 
laxes TT > 20 mas (i.e. G5-K3 dwarfs) and between (y^Sl 
and 1™2 for 3 < tt < 20 mas (i.e. G5-K2 giants and sub- 
giants). These criteria were imposed to select stars with a 
significant outer convective envelope where the likelihood 
of detecting magnetic activity is highest. Out of the 1,058 
stars, 371 (35%) were found with Call H&K emission but 
only 78 (7.3%) with v sin i > 10 km s^^. On the contrary, a 
lithium line was detected in 283 (74%) of all stars that had 
Call emission (with 58% of the stars with lithium above 
10 mA). Out of a subsample of 172 stars with moderate to 
strong Call emission, 168 (97.7%) turned out to be photo- 
metric variable and for 134 a photometric (rotational) pe- 
riod could be obtained. Finally, 36 targets were single-lined 
spectroscopic binaries (SBl), of which 17 were new detec- 
tions. A further 16 targets were found "possible SBls". An 
additional 30 targets were double-lined spectroscopic bina- 
ries (SB2), of which 19 were new detections. Two targets 
were even triple lined (SB3) of which one was a new de- 
tection and another four were new candidates. All along, 
there were a few misidentifications and misinterpretations 
as well as unrecognized literature entries. Whenever recog- 
nized, we try to clarify these in the present paper 

We now present time-series observations of all those tar- 
gets in the survey that were previously unknown or sus- 
pected spectroscopic binaries in 2000. This sample com- 
prised 59 stars presented in this paper. One additional star 
was added because of its comparable uncertain orbit despite 
being a fourth magnitude star (e UMi). During the final ob- 
servational stages for the present paper we learned that Grif- 
fin (^09\, Griffin ( |20T0l l and Griffin & Filiz Ak (|20T0| had 
picked up many of our original SB candidate stars from pa- 
per I that found their way into the new edition of the CABS 
catalog (Eker et al. |2008| l, and independently determined 
orbital elements. We compare our results whenever possi- 
ble. Stellar identifications and some basic observable prop- 
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Fig. 1 HD 171067 as an example for a STELLA/SES spectrum. The top panel shows the entire wavelength coverage from 
3880-8820 A. Note the increasing inter-order gaps starting around 7340 A. The lower panels zoom into three wavelength 
regions employed in this paper for analysis; from left to right, one of the five echelle orders that are used for the PARSES 
spectrum-synthesis fit (thick line), the Balmer Ha line, and the region around the neutral lithium line at 6708 A. HD 17 1067 
{V=l™2) is an inactive single, slightly evolved, G8 star with v sin i < 3 km s^^. Exposure time was 1200 s. 



erties for all our target stars are summarized in Table[T] Note 
that three of the stars actually do not show significant Ca II 
H&K emission and would not be dubbed magnetically ac- 
tive but are left in the sample because they were monitored 
initially in order to search for signs of binarity. All three 
turned out to be single stars though. Altogether, 18 of the 
SBl candidates were found to be single with some of them 
still members of a visual-binary system. These stars are also 
in Table [T] but are summarized in a separate table later in 
the paper In one case the two components of a wide vi- 
sual binary (BDh-1 1 2052AB = ADS 7406AB) were ti-eated 
as separate stars throughout the paper. One double-lined bi- 
nary (HD 16884) turned out to be actually a quadruple sys- 
tem with both pairs being SBl. First orbit determinations 
are presented for seven systems (HD 50255, HD 82841, 
HD 106855, HD 147866, HD 190642, HD 199967, and 
HD 226099). 

3 New observations and data reductions 

3.1 High-resolution optical spectroscopy 

Time-series high-resolution echelle spectroscopy was taken 
with the 1.2 m STELLA-I telescope between June 2006 and 
May 2012. Most spectra were exposed just long enough to 
measure a precise radial velocity and had S/N of between 
40-80: 1 but several spectra per target were exposed to reach 
S/N well above 100:1. A total of 6,609 spectra for a total of 
60 stars were obtained over the course of approximately six 



years. STELLA-I is a fully robotic telescope that, together 
with STELLA-II, makes up the STELLA observatory at the 
Izana ridge on Tenerife in the Canary islands (Strassmeier 
et al. 2004, 2010). The fiber-fed STELLA Echelle Spectro- 
graph (SES) is the telescope's only instrument. It is a white- 
pupil design with an R2 grating with two off-axis collima- 
tors, a prism cross disperser and a folded Schmidt camera 
with an E2V 2kx2k CCD as the detector. All spectra have 
a fixed format on the CCD and cover the wavelength range 
from 388-882 nm with increasing inter-order gaps near the 
red end starting at 734 nm towards 882 nm. The resolving 
power is i?=55,000 corresponding to a spectral resolution 
of 0.12 A at 650 nm (3-pixel sampling). An example spec- 
trum is shown in Fig. [T] We note that the SES received a 
major upgrade in summer 2012 with a new cross disperser, 
a new optical camera, and a new CCD. A bit earlier, the SES 
fiber was moved to the prime focus of the second STELLA 
telescope in 201 1. Further details of the performance of the 
system were reported by Weber et al. ( |2012 1 and Granzer et 
al. ( pOTO| ). 

SES spectra are automatically reduced using the IR^^ 
based STELLA data-reduction pipeline (Weber et al. 201 ] 
The images were corrected for bad pixels and cosmic-ray 
impacts. Bias levels were removed by subtracting the av- 
erage overscan from each image followed by the subtrac- 
tion of the mean of the (already overscan subtracted) master 



RA#} 
20T1) . 



' The Image Reduction and Analysis Facility is hosted by the Na- 
tional Optical Astronomy Observatories in Tucson, Arizona at URL 
iraf.noao.edu. 



www.an-joumal.org 



© 2012 WILEY- VCH Verlag GmbH & Co. KGaA, Weinhdm 



4 



K. G. Strassmeier, M. Weber, T. Granzer, S. Jarvinen: Rotation, activity, and lithium abundance in cool binary stars 



Table 1 Summary of program stars and name aliases. Listed are the Tycho V brightness and the B — V color in mag- 
nitudes, the spectral classification if available, the binarity according to this paper (S=single star, SB=spectrosopic bi- 
nary, SBl=single-lined SB, SB2=double-lined SB, SB3=triple-lined, VB=visual binary), the number of STELLA spectra 
^STELLA obtained and used for the orbit determination, its time span covered (At), and a primary reference if existent. 



HD 


Var. name 


HIP 


SAO 


V 


B — V 


Sp.type 


SB 


.^STELLA 


At 


Ref. 


HD553 


V741 Cas 


834 


11013 


8.17 


1.03 


KO 


SB2 


121 


1320 


s,34 




LN Peg 


999 


91772 


8.59 


0.81 


KO 


SB2+1 


135 


1302 


r 


HD8997 


EO Psc 


6917 


74742 


7.74 


0.96 


K2/K1-K6V 


SB2 


165 


1453 


b 


HD9902 


BG Psc 




74827 


8.71 


0.63 


F5-6V/G9-K0IV 


SB2 


95 


748 


a,c 


HD 16884 






110699 


8.94 


1.37 


K4III 


2xSBl 


83 


1046 


a,e 


HD 18645 


FUCet 


13968 


130230 


7.86 


0.75 


G2I1I-1V 


S 


70 


784 


f 


HD 18955 


IREri 


14157 


148731 


8.45 


0.82 


K0V/K2-3V 


SB2 


67 


1446 


g 


HD 23551 


MM Cam 


18012 


12924 


7.11 


0.91 


KOIII 


S 


101 


836 




HD 24053 




17936 


111446 


7.7 


0.8 


GO 


S 


77 


803 






AI Lep 




150676 


8.97 


0.57 


G6IV/G0V 


S 


93 


1460 


1 




HYCMa 




151224 


9.33 


1.00 


K0-3V-IV/K1V-IV 


SB2 


87 


1782 


1 


HD 40891 




28935 


13714 


8.40 


0.85 


G5 


SBl 


85 


1380 


h 


HD 43516 




29750 


151290 


7.37 


0.85 


G8I1I 


S 


57 


830 




HD 45762 


V723 Mon 


30891 


133321 


8.30 


0.87 


GO 


SB2+1 


89 


1231 


z 


HD 50255 




32971 


152024 


7.43 


0.68 


G4.5V 


SB2 


66 


1471 




HD 61994 




38018 


6310 


7.08 


0.67 


GO 


SB2 


105 


899 


t 


HD 62668 


BMLyn 


38003 


41995 


7.73 


1.10 


KOIII 


SBl 


137 


2000 


a,f,i 


HD 66553 




39515 


97536 


8.48 


0.85 


G5 


SBl 


92 


2001 


h 


HD 73512 




42418 


116990 


7.91 


0.90 


K0V/K4V 


SB2 


74 


892 


a 


HD 76799 




44007 


176747 


7.11 


0.99 


KOIII 


S 


75 


1442 








46634 


98614 


8.76 


0.86 


G5 


VB,S 


96 


1129 


y 


HD 82159 


GSLeo 


46637 


98615 


8.85 


0.92 


G9V 


VB,SB1 


107 


1458 


u,32 


HD 82286 


FFUMa 


46919 


14919 


7.89 


0.96 


G5 


SB2 


161 


1446 


i,31,35 


HD 82841 


OSHya 


46987 


136965 


8.45 


1.08 


K2 


SBl 


61 


1447 






EQ Leo 


50072 


99011 


9.39 


1.09 


Kim 


SBl 


62 


1169 


a 


HD 93915 




53051 


43492 


8.07 


0.68 


G5V/G6V 


SB2 


85 


1169 


a 


HD 95188 


XZLMi 


53747 


81610 


8.45 


0.74 


G8V 


S 


85 


891 


v,j 


HD 95559 


GZLeo 


53923 


81634 


8.83 


0.96 


K0V/K2V 


SB2 


113 


1162 


k,30 


HD 95724 


YYLMi 


54028 


62375 


8.96 


0.94 


G5V 


S 


76 


856 


V 


HD 104067 




58451 


180353 


7.93 


0.99 


K2V 


s 


76 


887 




HD 105575 


QYHya 


59259 


180519 


9.04 


0.93 


(K5/M1)/G4 


SB3 


51 


1164 


29 


HD 106855 


UVCrv 


59914 


180648 


9.59 


0.81 


KIV 


SB2,VB 


88 


1227 




HD 108564 




60853 




9.45 


0.98 


K5V 


S 


70 


867 




HD 109011 


NOUMa 


61100 


28414 


8.10 


0.94 


K2V 


SB2 


184 


1629 





HD 111487 


IM Vir 




138983 


9.69 


0.64 


G5 


SB2 


106 


1292 


p 


HD 112099 




62942 


119652 


8.23 


0.86 


KIV 


SBl 


80 


1142 


a 






63322 


63275 


9.27 


0.84 


G6V 


VB,SB1 


104 


937 


y,32 


HD 112859 


BQCVN 


63368 


44410 


8.09 


0.92 


F5V/K0in-IV 


SB2 


93 


1249 


a,w 




CDCVn 


63442 


44421 


9.39 


1.19 


KOIII 


SBl 


85 


1165 


a 


HD 120205 




67344 


158174 


8.3 


0.9 


G5 


S 


93 


849 




HD 127068 


HKBoo 


70826 


101044 


8.43 


0.89 


G8V/G5-8IV 


SB2 


99 


1257 


a 


HD 136655 




75132 


83780 


9.01 


0.9 


KO 


S 


122 


842 


z 


HD 138157 


OXSer 


75861 


101580 


7.14 


1.02 


KOIII 


SBl 


148 


1086 


u,32 




V381 Ser 


77210 


121177 


9.16 


0.83 


K2V 


SB2 


96 


1114 


33 


HD 142680 


V383 Ser 


77963 


101816 


8.71 


0.95 


K0-2V/K7V 


SB2 


90 


781 


a,v 


HD 143937 


V1044 Sco 


78708 


184077 


8.65 


0.91 


G9V/M0V 


SB2 


86 


1131 


1 


HD 147866 


V894 Her 


80302 


84343 


8.1 


1.1 


KO 


SBl 


115 


1269 




HD 150202 


GIDra 


81284 


30015 


7.97 


0.93 


KOIII 


SBl 


179 


1112 


a 


HD 153525 


V1089 Her 


83006 


46403 


7.88 


1.04 


KOV 


S 


196 


1349 


X 


HD 153751 


eUMi 


82080 


2770 


4.22 


0.88 


G5III 


SBl 


220 


2144 


d 


HD 155802 




84303 


141567 


8.51 


0.89 


K3V 


S 


114 


857 




HD 171067 




90864 


103819 


7.19 


0.696 


G8V 


S 


143 


857 




HD 184591 




96280 


104991 


7.37 


0.85 


G5 


S 


155 


842 


x,z 


HD 190642 


V4429 Sgr 


99042 


163264 


8.08 


0.99 


F6-8V/K1I1I-IV 


SBl 


196 


1274 




HD 197913 


OR Del 


102490 


106466 


7.55 


0.76 


G6V/G8V 


SB2 


194 


1848 


m 


HD 199967 




103454 




8.01 


0.58 


G5 


SB2 


121 


945 
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Table 1 (continued) 



HD 


Var. name 


HIP 


SAO 


V 


B-V 


Sp.type 


SB 


A'^STELLA 


At 


Ref. 


HD 202109 


CCyg 


104732 


71070 


3.20 


0.99 


G8III 


SBl 


150 


1288 


n 


HD 218739 




114385 


52754 


7.14 


0.61 


G5 


S 


107 


791 




HD 226099 




97640 


68946 


8.01 


0.783 


G5 


SB2 


209 


1334 




HD 237944 




53209 


27812 


9.36 


0.71 


(G8V+/G8V+)/? 


SB3 


119 


1431 


a,q 



References, a: Griffin (2009), b: Griffin (1987), c: Cutispoto et al. (2003), d: Climenhaga et al. (1951), e: Favata et al. (1993), f: Fekel 
(1997), g: Fekel et al. (2004), h: Latham et al. (2002), i: Henry et al. (1995), j: Wright et al. (2004), k: Karatas et al. (2004), 1: Cutispoto 
et al. (1999), m: Griffin (2005), n: Griffin & Keenan (1992), o: Halbwachs et al. (2003), p: Morales et al. (2009), q: Otero & Dubovsky 
(2004), r: Fekel et al. (1999), s: Duemmler et al. (2002), t: Duquennoy & Mayor (1988), u: Galvez et al. (2006), v: Barnes (2007), w: 
Montes et al. (2000), x: Strassmeier (1994), y: Galvez et al. (2005), z: Griffin (2010), 29: Szalai et al. (2007), 30: Galvez et al. (2009), 
31: Galvez et al. (2007), 32: Griffin & Fihz Ak (2010), 33: Goldberg et al. (2002), 34: Griffin (2003), 35: Griffin (2012). 



Table 2 STELLA radial velocities of standard stars. N is the total number of spectra, (wstella) is the average radial 
velocity from STELLA spectra, ctstella its standard deviation, and (f coravel) the CORAVEL velocity from Famaey et 
al. (|2005| or, if not in Famaey et al., from Udry et al. ( 1999[ l. The velocity difference STELLA-minus-CORAVEL is given 
in the last column. The averaged offsets from all stars weighted by the number of data points is given in the last line. All 
velocities are in kms^^. 

HD Name Class. (ustella) o-stella {«coravel) "stella-coravel 



HD 3712 


a Cas 


KO II-III 


230 


-3.863 


0.094 


-4.31 


+0.447 


HD4128 


/?Cet 


KOIII 


143 


+13.686 


0.107 


+13.1 


+0.586 


HD 12929 


a Ari 


K2III 


440 


-14.138 


0.118 


-14.64 


+0.502 


HD 18884 


a Cet 


M2 III 


60 


-25.688 


0.293 


-26.08 


+0.392 


HD 20902 


a Per 


F5Ib 


152 


-1.800 


0.138 


-2 


+0.200 


HD 25025 


7 Eri 


MO III 


26 


+61.159 


0.217 


+61.1 


+0.059 


HD 29139 


a Tau 


K5 III 


503 


+54.646 


0.188 


+54.26 


+0.386 


HD 36673 


a Lep 


FOIb 


53 


+25.871 


0.413 


+25.2 


+0.671 


HD 62509 


/?Gem 


KOIII 


442 


+3.783 


0.076 


+3.23 


+0.553 


HD 81797 


a Hya 


K3III 


203 


-4.033 


0.076 


-4.7 


+0.667 


HD 84441 


e Leo 


GO II 


165 


+5.001 


0.046 


+4.5 


+0.501 


HD 107328 


16 Vir 


K0.5 III 


332 


+37.117 


0.160 


+36.4 


+0.717 


HD 109379 


/3Crv 


G5 II 


49 


-6.947 


0.069 


-7.6 


+0.653 


HD 124897 


a Boo 


K1.5 III 


787 


-4.871 


0.149 


-5.3 


+0.429 


HD 146051 


5 0ph 


M0.5 Illab 


146 


-19.133 


0.146 


-19.6 


+0.467 


HD 156014" 


a Her 


M5Ib-II 


303 


-31.168 


0.694 


-32.09 


+0.922 


HD 161096 


/3 0ph 


K2 III 


763 


-11.844 


0.095 


-12.53 


+0.686 


HD 168454 


5Sgr 


K2III 


18 


-19.798 


0.035 


-20.4 


+0.602 


HD 186791 


7 Aql 


K3II 


538 


-2.430 


0.261 


-2.79 


+0.360 


HD 204867 


/? Aqr 


GO lb 


80 


+6.808 


0.151 


+6.3 


+0.508 


HD 206778 


ePeg 


K2Ib 


572 


+3.588 


0.262 


+3.39 


+0.198 


HD 212943 


35 Peg 


KOIII 


402 


+54.740 


0.111 


+54.16 


+0.580 



Weighted average +0.503 



"Double star with the A component being a "semiregular pulsating variable". 



bias frame. The target spectra were flattened by dividing by 
a nightly master flat which has been normalized to unity. 
The nightly master flat itself is constructed from around 
50 individual flats observed during dusk, dawn, and around 
midnight. After removal of the scattered light, the one- 
dimensional spectra were extracted with the standard IRAF 
optimal extraction routine (Horne |1986| l. The blaze func- 
tion was then removed from the target spectra, followed by 
a wavelength calibration using consecutively recorded Th- 
Ar spectra. Finally, the extracted spectral orders were con- 
tinuum normalized by dividing them with a flux-normalized 



synthetic spectrum of the same spectral classification as the 
target in question. 

3.2 VI and by photometry 

Johnson-Cousins V{I)c and/or Stromgren hy photometry 
of 56 of the targets were obtained with the Amadeus (T7) 
and/or the Wolfgang (T6) automatic photoelectric telescope 
(APT) at Fairborn Observatory in southern Arizona, respec- 
tively. These 56 stars are listed in Table |6] along with the 
analysis results. A total of 16,591 observations in either 
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V{I)c or by pairs were obtained. Typically, one APT obser- 
vation consists of three ten-second (30 s for Stromgren) in- 
tegrations on the variable, four integrations on the compari- 
son star, two integrations on the check star, and two integra- 
tions on the sky. A 30" diaphragm was used for all targets. 
The standard error of a nightly mean for Amadeus from 
the overall seasonal mean changed over the past decade but 
was typically 4-6 mmag in V and 6-8 mmag in Ic for 
the brightness range of stars in this paper. The observing 
seasons between 2007-2009 showed increasing scatter due 
to a slowly but systematic malfunction of the acquisition 
CCD camera. By late 2009 it got so bad that we had to ex- 
change the entire camera plus its CCD and after that, start- 
ing with HJD 2,455,143, the scatter was back to the original 
values quoted above. The standard error of a nightly mean 
for Wolfgang was between 1 .2-3 mmag, depending on sys- 
tem brightness. For further details we refer to Strassmeier et 
al. ^997} and Granzer et al. ( |200T] i. 

From concuiTent observations of Johnson standards in 
V{I)c, we also deduce an all-sky solution and apply it to 
the differential values whenever feasible. Its accuracy never 
significantly exceeds O'.'^Ol in V though. In case B data were 
taken, B — V can be determined only relative to the compar- 
ison star because our all-sky standards were not observed in 
the B band due to time constraints. Absolute errors are typi- 
cally around 0™0 1 for A{V — Ic) except for the time period 
2007-2009, as mentioned above, and was then likely 0'."02. 

4 Data products 

4.1 Radial velocity precision and zeropoint 

Radial velocities were determined from an order-by-order 
cross coiTelation with a synthetic template spectrum from an 
ATLAS9 atmosphere (Kurucz 1993) l that roughly matches 
the target spectral classification. A two-dimensional cross 
correlation is performed in case the target is double lined. 
Each of 19 selected echelle orders gives one relative radial 
velocity. These velocities are then weighted according to 
the spectral region depending on the amount of available 
spectral lines and are then averaged. The true internal er- 
ror is likely close to the rms of these relative velocities. 
Numerical simulations of the cross coiTelation errors from 
synthetic spectra and various two-Gaussian fits to its peaks 



are discussed in the paper by Weber & Strassmeier (201 1 1. 
Note that the external rms values were significantly larger 
during the initial year of STELLA operation in 2006/07 
(«120 ms^^), compared to thereafter (w30 ms^^). The 
final radial velocities of our program stars are barycentric 
and corrected for Earth rotation. No gravitational redshift 
corrections are applied. 

Table |2] lists our standard-star measurements. A total 
of 22 radial-velocity standards were monitored since mid 
2006. Their velocities were derived similarly from cross 
correlations with a synthetic template spectrum that fitted 
the star's classification. Among these standards are often 



Table 3 Average radial velocities for the single stars. At 
is the time span of STELLA observations in days from Ta- 
ble [T] is the average radial velocity and its standard 
deviation. Radial-velocity variations are detected for some 
of the stars due to stellar rotation and an asymmetrically 
spotted surface. Rotation periods, P, if existent, are given 
in days. 



Name 


At 




Note 




(days) 


(kms^^) 






/ 


o 1 ^ 1 _i_n n^"? 


P— 9 1 '>rl 


I iLJ ZjJJ 1 






i—o 1 . /a 


I I FA O/in^'^ 
I ILJ Z'+UJJ 








SAO 150676 


1460 


+25.69±0.53 


P=1.79d 


HD 43516 


830 


+21.846±0.059 




HD 76799 


1442 


-31.015±0.187 


P=50.13d 


HIP 46634 


1129 


+27.927±0.047 


Var. but no P 


HD 95188 


891 


+6.439±0.051 


P=7.00d 


HD 95724 


856 


+3.365±0.055 


P=11.49d 


HD 104067 


887 


+15.355±0.066 




HD 108564 


867 


+111.363±0.009 




HD 120205 


849 


-28.870±0.048 




HD 136655" 


842 


-31.921±0.035 




HD 153525 


1349 


-6.844±0.039 


P=15.4d 


HD 155802 


857 


-32.71 1±0.033 




HD 171067° 


857 


-45.771±0.042 




HD 184591° 


842 


-37.635±0.044 


P=48.9d:: 


HD 218739 


791 


-5. 201 ±0.040 





"Without Call H&K emission according to our paper I and thus 
not an active star. 



observed bright giant stars like j3 Oph, a Ari, j3 Gem, 
a Tau, 16 Vir, a.o., which are all lAU radial-velocity stan- 
dards (see Scarfe et al. I990| l. These spectra are also used 
to identify two epochs in our data when all stars on the ob- 
serving menu showed a constant offset in velocity. These 
epochs are identified with two of our hardware maintenance 
episodes where the fiber injection unit had been readjusted. 
The average offsets determined from above standard stars 
were —0.244 kms^^ and —0.366 kms^^ for the epochs 
2,453,902-4,085 and 2,454,192^35, respectively. These 
corrections were applied to all data in this paper. 

The preliminary zero point determination in Strassmeier 
et al. ( 20I0| l from four standard stars ifi Oph, a Cas, 
P Gem, /3 Cet) showed an offset of the STELLA spectra by 
H-0.46 km s~^ with respect to the CORAVEL system defined 
by Udry et al. ( |I999p . In the present paper, we revise this 
value to H-0.503 km s~^ from 22 standard stars. The weights 
applied for this average shift are the total number of spectra 
per radial- velocity standard. Note that an absolute heliocen- 
tric zero point for STELLA has not been determined yet, 
e.g. by using asteroids. All velocities in the present paper 
are on the STELLA system. 

Table [3] lists the stars in our sample that were found 
to exhibit no orbital-induced radial velocity variations, i.e. 
are presumably single stars. These stars are mostly still 
very active stars and further investigation with the present 
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P=1555.25 days 



HP 16884ab 




P=69.302 days 




0.4 0.6 
Phase 



P=3. 855871 days 




P=106.63 days 



0.4 0.6 
Phase 



0.4 0.6 
Phase 



0.4 



0.6 



Phase 
HP 105575(ab)c 



P=49.0823 days 




P=75.8996 days 



1/ 



P=56.732 days 




P=4550 days 




Fig. 2 Velocity curves for the STELLA single-lined orbital solutions. The dots are the observations and the lines are 
the orbital fits. Velocities are in kms^^ and the horizontal axes are in orbital phase. The horizontal dashed line in each 
plot indicates the systemic velocity. STELLA data are shown as filled dots. For HIP 999, the triangles are from Fekel et 
al. ( |1999| l. For HD 105575, values from paper I are shown as down-ward pointing triangles, the one FEROS point as a 
diamond, and the two HARPS points as upward-pointing triangles. For the other long-period system, HD 202109, data 



from Griffin & Keenan ( 1992 1 are added and shown as triangles. 
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epsUMi P=39.4804 days HP 190642 P=49.597 days 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Phase Phase 



HP 202109 p=6447.29 days 




0.0 0.2 0.4 0.6 0.8 1.0 

Phase 



Fig. 2 (continued). 
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SAO 151224 



P=4. 97769 days 
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/ % 
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P=59.936 days 



P=454.91 days 
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P=552.8 days 



P=128.241 days 
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— 40 









© ■ -OX. 

.0' 




/ 6 










q 

q 

q 











0.4 0.6 
Phase 



Fig. 3 Velocity curves for the STELLA double-lined orbital solutions. The primary velocities are plotted as filled dots, 
the secondaries as open circles. Otherwise as in Fig. [2] For HD 237944 the triangles mark the tertiary component. 
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Fig. 3 (continued). 
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Fig. 3 (continued). 



time-series data are forthcoming. Several of them, e.g. 
SAO 150676 or HD 76799, show periodic radial-velocity 
variations with amplitudes of the order of 1 kms^^, which 
we interpret to be due to stellar rotation. Few of our single- 
star survey targets in paper I were also revisited by Grif- 
fin ( [20T0| , most notably HD 136655 and HD 184591, for 
which we recommend to read his "Notes on the six constant- 
velocity stars" for their observational history. 

4.2 Binary orbits 

For single-lined (SBl) as well as double-lined (SB2) bina- 
ries we solve for the components using the general least- 



square fitting algorithm MP FIT (Markwardt 2009) 1. For so- 
lutions with non-zero eccentricity we use the prescription 



from Danby & Burkardt ( 1983[ ) to calculate the eccentric 
anomaly. Usually, we first determine the orbit for both com- 
ponents separately and then combine the two and give rms 
values for both. Eclipsing binaries were treated slightly sep- 
arately due to the Rossiter-McLaughlin effect in the radial- 
velocity curve. Some of the velocities around the two con- 
junctions were discarded by applying 3-0" clipping. Most 
of the orbits are sampled nearly perfectly due to robotic 



scheduling (Granzer et al. |2001| l but the unexpected high 
eccentricity in three systems left some phase gaps during 
periastron. The STELLA time coverage does not exceed 
5.9 years, typical values are around 3.3 years. Although 
partly compensated for by the high velocity precision and 
the dense sampling, systems with periods longer than, say, 
200 days, have comparably uncertain orbital periods. When- 
ever there were literature data for the eight systems with 
orbital periods in excess of wlOO d that expanded the to- 
tal time range and improved the orbit, these are included in 
the period determination. This was the case for HIP 999, 
HD 105575(ab)c and HD 202109. Otherwise, the orbits are 
computed from STELLA radial velocities alone. The result- 
ing elements and their errors are summarized in Table |4] and 
Table |5] for the single-lined and the double-lined systems, 
respectively. Assuming that the fit is of good quality, we 
derive the element uncertainties by scaling the formal one- 
sigma errors from the covariance matrix using the measured 
values. Tpor is a time of periastron or, if the orbit is cir- 
cular, a time of maximum positive radial velocity. Note that 
we give orbital periods as observed and not corrected for the 
rest frame of the system. 
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The computed radial-velocity curves are compared with 
our observed velocities in the series of panels in Fig. [2] 
and Fig. |3] for the single-lined orbits and the double-lined 
orbits, respectively. Note that two of the four triple sys- 
tems, HIP 999 and HD 105575, appear in both figures. 
The latter system is even triple lined. The third triple sys- 
tem, HD 237944, is also triple lined and consists of an 
SB2 in a short orbit but with the third component in a 
wide thousands-of-years orbit around the close pair The 
fourth triple system is HD 45762 where the secondary is 
again an SBl but with a period exactly 1/3 of the wide 
SB2 pair. HD 106855 is a close visual double with its pri- 
mary A-component to appear double lined, its bona-fide 
more massive Ab component being again an SBl for itself 
HD 16884 appears to be a quadruple system with two SBl 
pairs in a bound 106.65-d orbit. More comments to individ- 
ual systems are collected in Sect. |5] Note that these com- 
ments do not intend to give a comprehensive overview of 
CDS/Simbad-collected literature but are meant to draw at- 
tention to some specific system features. 

At the time of the start of our monitoring program in 
2006, a total of 17 stars were already either known binaries 
or had even an orbit computed, or at least had an orbital 
period known. However, we accumulated enough observa- 
tions for all systems to allow for orbital solutions just us- 
ing our own data for a certain epoch. This was motivated 
not only by consistency arguments but also by detections 
of orbital-period variations in active binaries and partic u 
larly in one of our targets (FF UMa; Galvez et al. |2007f 



is planned for the future whenever feasible (see, e.g. Strass- 



Orbital-period variations of magnetically active binaries are 
known for a long time (e.g. Hall & Kreiner [1973| see also 
Lanza et al. |2006 ) and were partly explained by mild mass 
exchange and mass loss due to a stellar wind. However, Ap- 
plegate ( 1992[ l proposed a connection between the gravita- 
tional quadrupole moment and the shape of a magnetically 
active component in the system as it goes through its mag- 
netic activity cycle. Therefore, it appears rewarding to rede- 
termine the orbits of such binaries from time to time. 



4.3 Rotational periods and spot amplitudes from 
photometry 

Table |6] is the summary of the results from our photomet- 
ric analysis. A total of 56 stars out of the total sample of 
60 targets were observed, of which 39 have a detectable 
period. The time coverage and the sampling are very dif- 
ferent from target to target, ranging from a minimum of 20 
nights for QY Hya (HD 105575) to 4,160 nights (>11 yrs) 
for HY CMa (SAO 151224) with a sampHng of all-night- 
long monitoring (QY Hya) to one observation per night as 
for most targets. This diversity of data wealth makes a co- 
herent spot-modelling analysis not feasible for the present 
paper and we just focus on the determination of average 
light-curve parameters. A more detailed analysis per target 



meieretal. 2011 for the most recent example ofHD 123351, 
excluded from this paper). 

Most relevant for this paper is the determination of a 
precise photometric period, Pphtm, that can be interpreted 
as the stellar rotation period. This allows a superior mea- 
surement of rotation with respect to spectroscopic v sin i 
measurements because they are independent of the (un- 
known) inclination of the rotational axis and can be deter- 
mined much more precisely. To proceed, we first pre-filtered 
all photometric data by excluding data points with an rms 
of >0.02 mag in order to remove data grossly affected by 
clouds. Periods were then searched for either in Johnson V 
or in Stromgren y or, for the stars where we have both pass- 
bands, in the combined y + V series. For light-curves with 
a highly harmonic content, the Lomb periodogram (Lomb 
|1976| l in the formulation of Scargle (1982) was applied. In 
cases where a period could have been affected by the win- 



See our individual notes for this star tliougli. 



dow function, the CLEAN algorithm (Roberts et al. 1987 1 
was employed to verify, but not to alter, the period identi- 
fication from the Lomb periodogram. For light curves with 
a highly non-harmonic content, i.e. for eclipsing binaries, 
a string-length minimization, a variant of a Lafler-Kinman 
(Lafler & Kinman |1965|l statistic described in Dworetsky 
( 1983| l, was used to determine the final photometric period. 
Some of the stars had already had a period determined in 
our paper I but the period in the present paper supersede 
these values. In Fig. [4] a phased light curve from all data for 
each of the variable targets is shown. These light curves are 
phased with the photometric period indicated in the graph 
and the respective time of periastron listed in Tables [4] or 
|5] or ascending node for circular orbits. Single-star zero 
epochs are arbitrary. 

Errors for the periods were obtained by using a method 
sometimes referred to as "refitting to synthetic data sets" 
(e.g. Ford 2005) . This method estimates confidence inter- 
vals by synthesizing a large number of data sets (typically 
10^) out from the original data by adding Gaussian random 
values to the measurements proportional to the actual rms 
of the data. The respective period-search algorithm is then 
applied to these synthetic data and the resulting standard de- 
viation assumed to be the standard deviation of the original 
period. This value is given in Table|6]as the error for Pphtm- 

Note that Table |6] strictly lists photometric periods. We 
interpret these usually as stellar rotational periods but em- 
phasize that care must be taken because, e.g., in case the 
star is a close binary, an ellipticity effect could mimic a spot 
light curve with a Poib/2 photometric period that is easily 
misinterpreted as the stellar rotation period. For one system, 
HD 138157, we interpret the 55-mmag amplitude partly due 
to the ellipticity effect. Its true rotational period is close to 
the orbital period. For five of our six eclipsing binaries the 
rotational modulation appears to have the same period as the 
orbital motion. For one case, HD 190642, a newly discov- 
ered eclipsing pair, the spot wave appears doubled humped 
and gives half of the true rotation period. 
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Table 4 STELLA single-lined orbital solutions. Column "rms" denotes the quality of the orbital fit for a point of unit 
weight in kms^^. 



Name Porb Tpcr K 7 e cj aisini f{M) rms 

(days) (HJD245+) (kms-^) (kms"^) (deg) (lO" km) 



HIP999(ab)c 


1555. 


0885.7 


5.99 


-0.10 







128. 


0.0348 


1.154 




±2.6 


±5.5 


±0.14 


±0.08 






±3.1 


1 rx rxrx n 

±0.0025 




HD 16884ab 


106.63 


4088.3 


29.97 


2.23 


0.0161 


266. 


43.93 


rx 

0.2978 


0.261 




±0.011 


±1.6 


±0.045 


±0.031 


±0.0014 


±5.2 


±0.066 


1 rx rxrx 1 ^ 

±0.0013 




HD 16884cd 


106.7 


4115.6 


16.3 


2.3 







24.0 


rx rx ^ 

0.0482 


2.215 




±0.17 


±0.77 


±0.38 


±0.26 






±0.56 


1 rx rxrxo a 

±0.0034 




HD 40891 


49.0823 


4092.901 


16.122 


17.791 


0.5105 


100.3 


9.357 


r\ rx 1 

0.01358 


0.080 




±0.00071 


±0.012 


±0.015 


±0.0098 


1 rx rxrxrx/^rx 

±0.00069 


±0.11 


1 rx rxrxrxo 

±0.0098 


1 rx r\r\r\r\ a 

±0.000043 




HD 62668 


69.3023 


4115.271 


31.10 


— 16.34 







29.64 


0.2165 


0.294 




±0.0016 


±0.022 


±0.037 


±0.028 






±0.036 


1 rx rxrxrxorx 

±0.00080 




HD 66553 


75.89957 


4110.790 


16.2 


34.564 


0.8719 


80.71 


8.30 


rx rxrxnrxr* 

0.00396 


0.130 




±0.00048 


±0.0075 


±0.19 


±0.015 


1 rx rxrxrxo 

±0.0028 


±0.32 


±0.13 


1 rx rxrxrx 1 

±0.00018 




HD 82159 


3.855871 


4064.9028 


30.68 


29.261 


0.2640 


215.2 


1.569 


rx rx 1 no 

0.01038 


0.242 




±0.000006 


±0.0025 


±0.036 


±0.025 


1 rx rxrx 1 r* 

±0.0012 


±0.25 


1 rx rxrx 1 rx 

±0.0019 


1 rx r\r\r\r\ a 

±0.00004 




T Tl~\ 0'^ /I 1 

HD 82841 


56.7316 


A i^r" A AO 

4064.48 


17.552 


A rx 1^0 r» 

40.782 


0.0874 


126.2 


13.640 


rx rxo 1 A r\ 

0.03149 


0.092 




±0.0017 


±0.11 


±0.017 


±0.013 


1 rx rxrx 1 rx 

±0.0010 


±0.72 


±0.013 


±0.00009 




H1P50072 


34.3157 


(~»r»rx 'V 

4220.7 


13.95 


11.97 


0.014 


301. 


6.582 


rx rxrxrxr"~7 

0.00967 


0.310 




±0.0027 


±1.7 


±0.060 


±0.042 


1 rx rxn a 1 

±0.0043 


±18 


±0.029 


1 rx rxrxrx 1 

±0.00013 




HD105575 


4550. 


4061. 


10.11 


20.03 


0.391 


186.8 


582. 


0.381 


0.160 




±52 


±13 


±0.09 


±0.075 


1 r\ rxrx f rx 

±0.0059 


±1.2 


±8.7 


±0.012 




HDl 12099 


23.50447 


4082.691 


18.017 


— 16.858 


0.4466 


21.90 


5.2105 


r\ rx 1 rx rx 

0.01023 


0.067 




±0.00013 


±0.0051 


±0.015 


±0.0079 


1 rx rxrxr\r'^ 

±0.00057 


±0.094 


1 rx rx rx ^ r* 

±0.0046 


±0.000027 




HIP63322 


4.593843 


4078.3032 


38.530 


—4.429 







rx J r\ 

2.4340 


r\ rxrx^rxrx 

0.02729 


0.119 




±0.000005 


±0.00056 


±0.017 


±0.012 






±0.0011 


±0.000036 




HIP63442 


38.7156 


4071.693 


15.56 


—45.538 







8.282 


0.01514 


0.266 




±0.0026 


±0.031 


±0.042 


±0.030 






±0.023 


1 rx rxrxrx 1 r\ 

±0.00012 




HD138157 


14.36891 


4150.70 


28.83 


-6.642 


0.0051 


83. 


5.696 


0.03574 


0.216 




±0.00010 


±0.40 


±0.047 


±0.018 


±0.0009 


±9.9 


±0.0093 


±0.00018 




HD147866 


114.050 


4234.12 


6.713 


-36.108 


0.5141 


100.5 


9.030 


0.002261 


0.076 




±0.0073 


±0.043 


±0.012 


±0.0075 


±0.0013 


±0.22 


±0.018 


±0.000014 




HDl 50202 


68.47160 


4162.339 


22.913 


-6.9620 







21.573 


0.08553 


0.083 




±0.00097 


±0.0081 


±0.0090 


±0.0064 






±0.0085 


±0.00010 




eUMi 


39.48042 


3933.392 


31.849 


-9.8623 







17.290 


0.1325 


0.121 




±0.00012 


±0.0047 


±0.013 


±0.0085 






±0.0068 


±0.00016 




HDl 90642 


49.5969 


4166.0 


40.20 


4.631 


0.0126 


314. 


27.42 


0.3346 


0.445 




±0.0012 


±0.73 


±0.045 


±0.033 


±0.0011 


±5.3 


±0.031 


±0.0011 




HD202109'= 


6446. 


2421363. 


3.30 


16.41 


0.244 


40.7 


284. 


0.0220 


0.413 




±14 


±82 


±0.057 


±0.043 


±0.016 


±3.3 


±5.1 


±0.0012 





"Time of periastron, or ascending node for circular orbits. 

''Triple system. Orbit (ab) around c is given. 

'^Quadruple system with two SB 1 pairs in a synchronized orbit. 

''SBl orbit of the tertiary component, c, around the close (eclipsing) ab pair 

■^SBl with long period. The orbit given is from a combination of STELLA data and selected published values. Tpcr is 2421363. See 
individual notes. 



The typical error of the magnitude zero point is around 
0.01-0.02 mag, but never higher than 0.05 mag. The wave 
amplitude, AV, in Table|6]is the maximum observed ampli- 
tude within our data coverage. For the eclipsing systems, we 
cut out the data points during eclipse phases and then rede- 
termined the maximum amplitude. The color index, C.L, is 
the average color index over all observations. In case of T7 
observations it is (V-I)c, for T6 it is b-y, if available. Vari- 
ations in the C.I. are present but are in general too noisy for 



period analysis. Most of the observations of HD 136655 had 
to be discarded due to obvious instrumental problems. Its 
remaining data shows no evidence for photometric variabil- 
ity. Orbital period and Tq of the eclipsing binary HD 105575 
are equally constrainable from photometry than from radial 
velocities. If e=0 is adopted, we obtain ro=2,454,160.772 
and Porb=0.29238±0.00009 d just from our 20-day long 
time series. However, because the photometry was obtained 
around JD2,451,250, i.e. almost 10,000 eclipses prior to the 
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Table 5 STELLA double-lined orbital solutions. Suffices 1 and 2 denote the primary and the secondary star, respectively. 



Name 


Period 


^ Per 


Ki 


T 
/ 


g 




(i\ sin i 


Ml sin^ i 


rmsi 








K2 








C12 sin i 


M2 sin"^ i 


rms2 






(HJD245+) 


(kms~^) 


(kms~^) 






(10® km) 




(kms~^) 


HD553 


9.059792 


4087.1233 


63.312 


-38.014 







7.8874 


1.019 


0.411 




±0.000032 


±0.0019 


±0.037 


±0.027 






±0.0047 


±0.0073 










65.46 








8.155 


0.9854 


0.390 








±0.23 








±0.029 


±0.0037 




HIP 999'' 


1.8441405 


4063.6374 


62.40 


-14.24 







1.582 


0.447 


0.519 




±1.610"^ 


±0.00055 


±0.063 


±0.041 






±0.0016 


±0.0071 










94.77 








2.403 


0.2945 


3.292 








±0.68 








±0.017 


±0.0026 




HD8997 


10.98336 


4072.091 


39.124 


21.546 


0.0412 


195.2 


5.9040 


0.3956 


0.365 




±0.00004 


±0.029 


±0.027 


±0.021 


±0.0007 


±0.93 


±0.0041 


±0.0011 










47.00 








7.092 


0.3294 


0.400 








±0.060 








±0.0090 


±0.00063 




HD9902 


25.3722 


4079.516 


29.46 


56.956 


0.5068 


74.33 


8.861 


0.1813 


0.201 




±0.00026 


±0.0060 


±0.035 


±0.014 


±0.00074 


±0.11 


±0.012 


±0.00066 










30.22 








9.089 


0.1767 


0.154 








±0.046 








±0.015 


±0.00057 




HD 18955 


43.32161 


4086.710 


53.83 


31.273 


0.75874 


174.84 


20.89 


0.9982 


0.237 




±0.00022 


±0.0045 


±0.090 


±0.024 


±0.00050 


±0.080 


±0.039 


±0.0052 










61.02 








23.68 


0.8807 


0.480 








±0.12 








±0.053 


±0.0042 




SA0151224 


4.977691 


4069.9382 


68.07 


42.67 







4.659 


1.130 


1.130 




±0.000012 


±0.0021 


±0.085 


±0.057 






±0.0058 


±0.0041 










88.92 








6.087 


0.8652 


0.869 








±0.14 








±0.0099 


±0.0026 




HD 45762'= 


59.9363 


4110.40 


65.447 


2.23 


0.0150 


84.6 


53.935 


0.599 


0.533 




±0.0016 


±0.61 


±0.064 


±0.048 


±0.0009 


±3.7 


±0.053 


±0.058 










14.9 








12.3 


2.63 


6.4 








±1.1 








±0.87 


±0.069 




HD 50255 


454.910 


4152.91 


11.77 


19.059 


0.3440 


324.42 


69.11 


0.5414 


0.150 




±0.038 


±0.28 


±0.023 


±0.011 


±0.0011 


±0.23 


±0.14 


±0.0045 










16.89 








99.23 


0.3771 


0.176 








±0.063 








±0.37 


±0.0022 




HD 61994 


552.82 


4397.45 


10.46 


-22.038 


0.4230 


222.3 


72.07 


0.465 


0.118 




±0.25 


±0.26 


±0.018 


±0.0088 


±0.0012 


±0.18 


±0.013 


±0.0044 










15.81 








108.9 


0.3120 


0.117 








±0.065 








±0.46 


±0.0019 




HD 73512 


128.241 


4145.44 


23.200 


33.478 


0.26240 


277.27 


39.478 


0.7722 


0.087 




±0.0046 


±0.045 


±0.014 


±0.0086 


±0.00045 


±0.13 


±0.025 


±0.0016 










26.349 








44.836 


0.6799 


0.170 








±0.025 








±0.043 


±0.0011 




HD 82286 


3.2751924 


4066.0088 


29.775 


-1.473 







1.3410 


0.19296 


0.025 




±4 10-'' 


±0.00012 


±0.0043 


±0.0026 






±0.00019 


±0.00010 










64.282 








2.8951 


0.089378 


0.174 








±0.014 








±0.00062 


±0.00003 




HD 93915 


222.965 


4182.78 


22.35 


-17.288 


0.3795 


37.80 


63.40 


0.8732 


0.180 




±0.022 


±0.086 


±0.033 


±0.015 


±0.00088 


±0.16 


±0.097 


±0.0032 










23.10 








65.54 


0.8447 


0.196 








±0.037 








±0.11 


±0.0030 








4U/0.8iUU 


iUo.ooa 


4.ziz 


U 




Z.z ( 4Uz 


U.ozoz 


U.lUo 




±7.510-* 


±0.00003 


±0.011 


±0.0064 






±0.00022 


±0.0004 










109.93 








2.3068 


0.81644 


0.165 








±0.025 








±0.0005 


±0.00025 




HD105575'^ 


0.2923378 


4160.768 


87.4 


27.4 







0.351 


0.288 


8.6 




±1.710"® 


±0.0044 


±4.5 


±2.3 






±0.018 


±0.027 










158. 








0.635 


0.16 


13 








±6 








±0.024 


±0.016 
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Table 5 (continued) 



Name 


Period 


rpa 

Per 


Ki 


1 


e 




ai sin i 


Ml sin^ i 


rmsi 








K2 








a2 sin i 


M2 sin^ i 


rms2 




(days) 


(HJD245+) 


(kms^^) 


(kms^^) 




(deg) 


{10^ km) 


(M0) 


(kms^^) 


HD 106855 


0.6709843 


4097.2944 


41.8 


-1.1 







0.385 


0.0166 


1.811 




±0.0000026 


±0.0031 


±0.65 


±0.42 






±0.006 


±0.0012 










37.7 








0.348 


0.0183 


9.92 








±1.4 








±0.013 


±0.0009 




HD109011 


1274.7 


5213.63 


9.582 


-9.8616 


0.5071 


247.36 


144.8 


0.5130 


0.094 




±0.86 


±0.23 


±0.014 


±0.0044 


±0.00060 


±0.13 


±0.24 


±0.0021 










12.47 








188.4 


0.3941 


0.147 








±0.021 








±0.35 


±0.0014 




HDl 11487 


1.3086147 


4072.3403 


92.350 


13.17 







1.6618 


0.9217 


0.268 




±0.0000002 


±0.00014 


±0.027 


±0.021 






±0.0005 


±0.0028 










133.4 








2.401 


0.6381 


1.776 








±0.19 








±0.0033 


±0.0011 




HD 112859 


18.49896 


4072.681 


43.466 


22.267 







11.057 


0.9293 


0.300 




±0.00014 


±0.0038 


±0.031 


±0.020 






±0.0078 


±0.0027 










52.57 








13.37 


0.7683 


0.265 








±0.07 








±0.018 


±0.0015 




HDl 27068 


15.15011 


4108.69 


22.09 


-50.188 


0.0098 


161. 


4.601 


0.1377 


0.190 




±0.00011 


±0.23 


±0.026 


±0.014 


±0.00087 


±5.5 


±0.0055 


±0.00069 










31.06 








6.471 


0.09790 


0.364 








±0.071 








±0.015 


±0.00034 




HIP77210 


9.943224 


4106.9724 


42.096 


1.609 


0.06069 


83.32 


5.7451 


0.4962 


0.052 




±0.000011 


±0.0041 


±0.007 


±0.0047 


±0.00015 


±0.15 


±0.001 


±0.00060 










53.262 








7.2691 


0.39213 


0.203 








±0.030 








±0.0041 


±0.00028 




HD142680 


24.53381 


4169.6573 


37.257 


-83.830 


0.3158 


324.81 


11.926 


0.7375 


0.042 




±0.00006 


±0.0022 


±0.006 


±0.0041 


±0.00015 


±0.031 


±0.0021 


±0.0010 










47.398 








15.172 


0.5797 


0.205 








±0.031 








±0.010 


±0.00047 




HD143937 


0.91483076 


4150.5197 


118.56 


-55.307 







1.4914 


0.8036 


0.285 




±0.0000002 


±0.0001 


±0.04 


±0.030 






±0.0005 


±0.0019 










133.5 








1.679 


0.7137 


1.258 








±0.16 








±0.0020 


±0.0010 




HD197913 


1778.5 


6285.5 


10.941 


-26.969 


0.2600 


342.57 


258.4 


0.9649 


0.100 




±1.3 


±2.0 


±0.010 


±0.0060 


±0.00072 


±0.24 


±0.30 


±0.0026 










11.52 








272.1 


0.9161 


0.161 








±0.013 








±0.37 


±0.0021 




HDl 99967 


6.1002208 


4169.06 


65.12 


-6.830 


0.00244 


170. 


5.463 


1.177 


0.102 




±0.000005 


±0.06 


±0.076 


±0.007 


±0.00015 


±3.5 


±0.0064 


±0.0032 










83.89 








7.037 


0.9139 


0.149 








-1-0 1 n 








-1-0 nns^i 

ztzU.UUoO 


-l-fl 009/1 




HD226099 


18.78209 


4162.0653 


38.528 


7.5772 


0.30961 


242.63 


9.4618 


0.47739 


0.086 




±0.000019 


±0.0016 


±0.0068 


±0.0042 


±0.00014 


±0.031 


±0.0017 


±0.00032 










42.968 








10.552 


0.42807 


0.182 








±0.013 








±0.0033 


±0.00021 




HD237944'= 


5.5076262 


4070.103 


76.237 


-15.25 


0.0141 


105. 


5.7733 


1.042 


0.264 




±0.0000047 


±0.019 


±0.048 


±0.017 


±0.0003 


±1.3 


±0.0036 


±0.0015 










77.410 








5.8621 


1.027 


0.223 








±0.048 








±0.0036 


±0.0014 





"Time of periastron, or ascending node for circular orbits. 

SB2, but with a third component determined from the 7-residuals. The third-component orbit is given in the SBl list in Tabled See 
also the individual notes. 

'^Hierarchial triple system with lines from components a and c. The orbit ab is SBl with P=19.933 d. See individual notes. 

SB 3 with a third component in the spectrum denoted c. Its orbit around the (eclipsing) ah pair is given in Tabled 
^SB3. The third component, c, has a median velocity of -12.25±1.3(stdev) km s~^ for the time of our observations. 
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Table 6 Summary of photometric data and results. Listed are the comparison ("cmp") star used for the differential 
observations, the start date of observations ("HJD start"), the time span of observations (At) in days, the telescope and 
photometric bandpasses used (Johnson-Cousins VI for the T7 APT and Stromgren by for the T6 APT), and the number of 
actual data points N, each the mean of three individual readings. The results are summarizes in the subsequent columns 
and list Vmax, the maximum, i.e. least spotted, brightness during our observations; AV, the largest observed amplitude due 
to spots and, if eclipsing, due to primary eclipse; (C.I.), the average color index (V-I)c for T7 observations and/or (b-y) for 
T6; ^phtm> the photometric period and its likely error obtained from a refitting to synthetic data; and a note if applicable. 



HD/name 


HD/cmp 
star 


HJD start 
(2,4+) 


At 


Tel. 


TV 


Vmax 

(mag) 


AVmax 

(mag) 


<C.I.> 
(mag) 


(days) 


Note 


HD553 


443 


51107 


253 


T6+7 


217 


8.15 


0.10 




Wforb 


rotation 


HD553 














0.28 




9.061±0.002 


eclipsing 


HIP 999 


977 


50395 


700 


T6 


134 


8.53 


0.08 


0.56 


1.84475±0.00008 




HD 8997 


9472 


51098 


95 


T6 


69 


7.75 


0.02 




10.54±0.11 




HD 9902 


H1P7910 


51098 


112 


T6 


122 


8.68 


0.05 




7.445±0.02 




HD 16884 


17030 


51434 


3327 


T6+7 


286 


8.87 


0.16 


1.58 


65.44±0.05 




HD 18645 


18668 


51098 


727 


T6 


241 


7.85 


0.04 


0.55 


21.5±0.1 




HD 18955 


19223 


51099 


50 


T6 


29 


8.46 








var at 2.3cr 


HD 23551 


24513 


51447 


3729 


T6+7 


242 


7.08 


0.035 


0.94 


81.69±0.03 


P from T6 


HD 24053 


23570 


55570 


74 


T7 


130 


7.72 


0.06 


0.55 






SAO 150676 


37792 


51434 


745 


T6+7 


389 


9.09 


0.08 


0.79 


1.7889±0.0003 




HD 40891 


42756 


55570 


60 


T7 


71 


8.42 


0.06 


0.76 






SAO 151224 


43306 


51079 


4160 


T7 


1352 


9.27 


0.39 


1.27 


~ -^orb 


rotation 


SAO 151224 














0.68 




4.98242±0.00003 


eclipsing 


HD 43516 


44396 


55576 


67 


T7 


82 


7.40 


<0.03 


0.85 






HD 45762 


45168 


55576 


49 


T7 


61 


8.35 


0.15 


0.98 


62.0±1.6 


+ellipticity 


HD 50255 


50674 


55576 


49 


T7 


38 


7.48 


0.05 


0.62 






HD 61994 


63748 


55576 


49 


T7 


43 


7.05 


0.03 


0.64 






HD 62668 


64106 


51457 


3782 


T7 


568 


7.46 


0.40 


1.18 


67.470±0.007 




HD 66553 


65257 


55576 


29 


T7 


6 


8.52 


<0.02 


0.72 






HD 73512 


73764 


55572 


92 


T7 


131 


7.92 


<0.02 


0.80 






HD 76799 


77361 


55572 


92 


T7 


133 


7.25 


0.06 


0.92 






HD 82159 


82410 


51242 


89 


T6 


52 


7.97 


0.05 




3.055±0.003 


polluted 


HD 82286 


82719 


51459 


3780 


T7 


609 


7.84 


0.16 


1.12 


3.27629±0.00003 




HD 82841 


83046 


51476 


3763 


T7 


914 


8.22 


0.07 


1.15 


56.58±0.04 




HIP 50072 


88853 


51242 


96 


T6 


58 


9.44 


0.12 




33.782±0.009 




HD 93915 


95379 


55572 


138 


T7 


156 


8.10 


0.20 


0.65 






HD 95188 


94996 


51249 


102 


T6 


77 


8.51 


0.025 




7.00±0.06 




HD 95559 


95242 


51551 


3686 


T6+7 


1061 


8.83 


0.08 


0.93 


1.517000±410-® 




HD 95724 


96778 


51239 


112 


T6 


78 


8.97 


0.05 




11.49±0.02 




HD 104067 


104414 


55572 


168 


T7 


158 


7.98 


0.05 


0.92 






HD 105575 


105814 


51239 


20 


T6 


101 


8.98 




0.62 


~ -Porb 


rotation 


HD 105575 














0.36 




0.29238±0.00009 


eclipsing 


HD 106855 


106991 


51551 


3685 


T7 


992 


9.35 


0.11 


1.27 


0.6705258±6 lO"'' 




HD 108564 


108599 


55572 


168 


T7 


203 


9.50 


0.06 


1.02 






HD 109011 


109894 


51225 


135 


T6 


60 


8.10 


0.02 




8.4±0.2 




T-fn 1 1 1 487 


1 1 1 976 


S1 14.9 


41 1 f\ 


T7 
1 / 


1 S71 


Q 69 




\J.O\J 


^ J orb 


lUlclllUIl 


HD 111487 














0.76 




1.30859±0.000007 


ecUpsing 


HIP 63322 


113168 


51240 


65 


T6 


35 


8.89 


0.08 




2.227±0.003 


FAP 20% 


HD 112099 


111733 


55569 


171 


T7 


171 


8.22 


0.04 


0.70 






HD 112859 


112220 


51225 


135 


T6 


76 


8.23 


0.06 




18.34±0.05 




HIP 63442 


113828 


51240 


120 


T6 


80 


9.39 


0.12 




40.29±0.08 




HD 120205 


120969 


55572 


168 


T7 


199 


8.35 


0.06 


0.75 


14.3±0.7 




HD 127068 


126583 


51537 


3700 


T7 


972 


8.47 


0.18 


1.05 


18.187±0.0001 




HD 136655 


136643 


50492 


355 


T7 


465 


8.81 


<0.02 


1.00 




constant 


HD 138157 


138085 


51540 


3697 


T6+7 


1159 


7.11 


0.14 


1.10 


7.18357±0.00005 


P/2+ell. 


HIP 77210 


140750 


51260 


100 


T6 


49 


9.14 


0.02 




13.7±0.1 




HD 142680 


142245 


51260 


100 


T6 


64 


8.70 


0.03 




33.4±1.5 




HD 143937 


143438 


51559 


48 


T7 


57 


8.67 


0.10 


1.16 




rotation 


HD 143937 














0.41 




0.91479±0.00003 


eclipsing 


HD 147866 


147266 


51447 


728 


T6+7 


92 


8.08 


0.03 


1.14 


80.3±2.2 
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Table 6 (continued) 



HD/name 


HD/cmp 


HJD start 


At 


Tel. 


N 




ATT" 

AVmax 


<C.I.> 


^phtm 


Note 




star 


(2,4+) 








(mag) 


(mag) 


(mag) 


(days) 




HD 150202 


149843 


55572 


168 


T7 


210 


7.94 


0.05 


0.73 


36.4±0.5 


1/2 Prot 


HD 153525 


153286 


51260 


100 


T6 


53 


7.94 


0.02 




15.4±0.1 




HD 155802 


156227 


55580 


160 


T7 


149 


8.50 


0.03 


0.85 






HD 171067 


170651 


55602 


138 


T7 


194 


7.17 


0.02 


0.55 






HD 184591 


183849 


51099 




T6 


263 


7.34 


<0.01 


0.62 


(48.9±2.2) 


uncertain 


HD 226099 


188149 


55634 


106 


T7 


210 


7.98 


0.05 


0.65 






HD 190642 


189893 


51434 


3668 


T7 


791 


8.15 


0.32 


1.26 




rotation 


HD 190642 














0.43 




24.8339±0.0009 


eclipsing 


HD 197913 


198109 


51098 


51 


T6 


33 


7.13 


0.02 




6.556±0.07 




HD 199967 


199354 


55652 


88 


T7 


160 


7.65 


0.03 


0.40 






HD 237944 


93470 


55569 


169 


T7 


675 


9.27 


0.07 


0.50 


5.41±0.013 


rotation 


HD 237944 














0.25 




5.5075±0.0003 


eclipsing 



spectroscopic data, the combination of the two data sets still 
results in a loss of the cycle count and could not be com- 
bined. Table |5]gives the orbit computed only with the radial 
velocities. 

4.4 Rotational velocities 

Selected spectra of our program stars were subjected to a de- 
noising procedure based on a principal component analysis 



developed for Doppler imaging (Carroll et al. 2007 2009 1. 



This procedure employs typically a total of around 1100 
spectral lines in the wavelength range 480-850 nm to deter- 
mine the noise spectrum. This noise spectrum can then be 
used to de-noise an arbitrary wavelength section. The only 
free parameter is the number of principal components. No 
strict rule can be given for the choice of this value but we 
followed the prescription laid out in Martinez Gonzalez et 
al. ( 2008| l. It suggests 15 components for our typical late 
spectral types and wavelength coverage. It is then used to 
de-noise a well-known and unblended spectral line from 
which we measure wsini under the assumption of a tem- 
perature and gravity dependent macroturbulence. For this 
paper, we have chosen the well-known and unblended Fe I 
549.7516 nm line with an excitation energy of 1.011 eV, a 
transition probability of log gf = —2.849, and a typical mi- 
croturbulence of 2.0 kms~^. The radial-tangential macro- 
turbulences, Crt, are taken from Gray (2005), as listed in 
Fekel ([1997]). This macroturbulence is then subtracted from 
the measured total line broadening. At this point, we note 
that Gray (|2005)l gave most probable macroturbulent veloci- 
ties, which are 1 .414 times greater than the rms values given 
here. The total line broadening is measured independently 
from fits with synthetic spectra over a large range of wave- 
lengths (see next section) and the resulting v sin i measures 
were then compared and found to agree within their esti- 
mated errors. The final values of v sini and (^rt for all our 
program stars are tabulated in Table |7]for further reference. 

Up to three nightly (and two daily) Th-Ar calibration 
frames are used to monitor the spectrograph focus through- 
out the year. The line widths of Th-Ar emission lines are 



determined automatically and an unweighted average stored 
in the STELLA data base. No focus drifts were evident dur- 
ing the time in question. However, manual refocusing was 
done after every maintenance run. The telescope focus is be- 
ing controlled and adjusted daily by insetting a focus pyra- 
mid into the telescope beam. It splits a stellar image into 
four equally distant images if the telescope is in focus (see 
Granzer et al. |2010| l. The STELLA control system issues a 
message to the operator in Potsdam if this had to be cor- 
rected. 

4.5 Stellar atmospheric parameters 

Five selected spectral orders (#28, 29, 33, 37, 39) from 
the STELLA/SES spectra are used to determine the stel- 
lar effective temperature, the gravity and the metallicity. 
These orders cover the wavelength ranges 549-556, 560- 
567, 583-591, 608-616, and 615-623 nm. Our numeri- 
cal tool PARSES ("PARameters from SES"; AUende-Prieto 
[2004') is implemented as a suite of Fortran programs in the 
STELLA data analysis pipeline. It is based on the synthetic 
spectrum fitting procedure described in Allende-Prieto et 
al. ('2006'). Model atmospheres and synthetic spectra were 
taken from the ATLAS9 CD (KuiTiczfl993])- Synthetic spec- 
tra are pre-tabulated with metallicities between —2.5 dex 
and +0.5 dex in steps of 0.5 dex, logarithmic gravities be- 
tween 1.5 and 5 in steps of 0.5, and temperatures between 
3500 K and 7,250 K in steps of 250 K for a wavelength 
range of 380-920 nm. All calculations were done with a 
microturbulence of 2 km s~ ^ . This grid is then used to com- 
pare with the five selected echelle orders of each spectrum. 
Free parameters are Tcff, log 5, [Fe/H], and line broaden- 
ing ^yivsini)^ + Cprp, where Crt is the radial-tangential 
macroturbulence, again adopted from Gray (2005 ) and was 
fixed during the fit. Internal errors are estimated from the 
rms of the solutions to the five echelle orders and were for 
Tcs typically 20-30 K, for log g typically 0.06 dex, and for 
[Fe/H] typically 0.03 dex. 

Double-lined spectra can not be run automatically 
through the PARSES routines. We first apply a similar dis- 
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Fig. 4 Light curves for the stars in Table|6] The time coverage of the observations. At, and the photometric period, P, are 
indicated on the top of each panel, both in days. The dots are the F or y observations from the full epochs given in Table|6] 
Error bars indicate the standard deviation from three individual measurements. Photometric phase has been computed with 
the photometric period and the respective time of periastron or ascending node for circular orbits. Single star light curves 
are plotted with an arbitrary starting time. Note that most of the scatter is due to intrinsic spot changes and not due to 
instrumental origin. The standard error of a mean was 4-6 mmag for the V data and 1.2-3 mmag for the y data. 
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Fig. 4 (continued) 



entangling technique as for the automatic radial-velocity 
measurement of triple-lined spectra, i.e. we remove once 
the primary and once the secondary component from the 
combined spectrum at well-separated phases in order to cre- 
ate a single-star spectrum for each component. The result- 
ing spectra were then subjected to the PARSES synthetic- 
spectrum fitting. This was done only for the five echelle or- 
ders that are employed for the PARSES fit (see above) and 
not for the entire 80 orders. Errors are generally larger for 
SB2s than for SBl and single systems due to continuum 
cross talk. External errors are again estimated from the rms 
of the solutions to the five echelle orders and were for Tas 
typically 100 K, for log 5 typically 0.3 dex, and for [Fe/H] 
typically 0.3 dex. However, the range of errors is systemati- 
cally larger for the SB2 systems than for the SB 1 and single 
stars because of the large range of brightness differences 
between the components in our sample (ranging from 
i.e. nearly equal components, to 12:1 for HIP 999). 

Lithium abundances are determined for all individual 
spectra that had sufficient S/N in the respective echelle or- 



der. The time averaged results are listed in Table |7] with re- 
spect to the log n(H) = 12.00 scale for hydrogen. A double 
Gaussian fit to Li I 6707.8 A and the nearby Fel 6707.4- 
A line yields an individual equivalent width for the lithium 
line. Despite that these values are still the combined equiv- 
alent width for the close blend from the ^Li and the '^Li 
isotopes, it effectively removes the Fe I blend. Weaker CN 
features are neglected in this study because of limited S/N 
and spectral resolution. The listed equivalent widths are the 
averages of the individual measures with unit weight. The 
errors in Table [T] are still the O-C error from an individual 
Gaussian fit but again averaged over all measures in time. 
Because the rms values of the equivalent widths of an entire 
time series per star were never larger than 3 ct of above O- 
C error, we do not explicitly list them in Table [7] However, 
these rms are useful to judge eventual long-term variabil- 
ity. SB2's are measured only at phases of quadrature, where 
blending is minimized. The curves of growth from Pavlenko 
& Magazzu ( |1996| l are then used to convert the equivalent 
width to a logarithmic lithium abundance using the non-LTE 
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Table 7 Determination of astrophysical properties. Column SB denotes the binary component in case the system is 
SB2 or even SB3 (S means single star, a means primary component, b secondary component, c tertiary component). The 
parallaxes are from the revised Hipparcos catalog (van Leeuwen |2007| l. Effective temperatures, gravities and metallicities 
are from the STELLA PARSES analysis. Note that errors are mean internal errors from the fits to five spectral orders. The 
Li abundance is from equivalent-width measures of the Li I 670. 8nm line converted with the NLTE calibration of Pavlenko 
& Magazzu ( |1996 l and the temperatures Usted here. Crt is the radial-tangential macroturbulence adopted. 



Name SB tt T^a logp [Fe/H] logn(Li) ("rt tjsini 

(mas) (K) (ms"2) (mA) (H = 12.00) (kms^^) 



HD 18645 


S 


8.71 


5445±30 


3.53±0.04 


-0.19±0.03 


14±4 


1.47±0.15 


3 


11.0±0.2 


HD 23551 


S 


5.00 


5045±25 


3.04±0.05 


-0.11 ±0.02 


12±4 


1.03±0.16 


3 


6.4±0.2 


HD 24053 


S 


30.74 


5600±25 


4.09±0.03 


-0.01 ±0.02 


49±8 


2.26±0.12 


3 


<3 


SAO 150676 


S 




5750±70 


4.54±0.11 


-0.16±0.07 


97±10 


2.71±0.10 


3 


23.6±0.7 


HD 43516 


s 


4.10 


5120±50 


2.86±0.07 


-0.15±0.04 


39±6 


1.72±0.15 


3 


<3 


HD 76799 


s 


6.13 


4915±25 


3.26±0.07 


-0.12±0.03 


12±4 


0.86±0.18 


3 


7.8±0.3 


HIP46634 


s 


27.46 


5140±35 


4.44±0.04 


-0.01 ±0.03 


33±6 


1.59±0.14 


2 


<3 


HD 95188 


s 


27.63 


5330±40 


4.39±0.06 


-0.18±0.03 


13±6 


1.29±0.24 


3 


6.9±0.7 


HD 95724 


s 


28.37 


4960±20 


4.45±0.06 


+0.06±0.02 


15±5 


0.97±0.16 


2 


<3 


HD 104067 


s 


48.04 


4820±20 


4.37±0.04 


-0.03±0.02 


16±4 


0.83±0.15 


2 


8.1±0.3 


HD 108564 


s 


35.30 


4560±25 


4.40±0.06 


-0.90±0.03 


21±7 


0.66±0.20 


2 


14.8±1.4 


HD 120205 


s 


31.78 


5260±30 


4.36±0.05 


+0.03±0.03 


10±5 


1.13±0.18 


2 


<3 


HD 136655 


s 


24.87 


5010±25 


4.27±0.03 


+0.09±0.04 


11±5 


0.92±0.20 


3 


<3 


HD 153525 


s 


57.13 


4775±25 


4.47±0.05 


-0.12±0.03 


12±5 


0.66±0.20 


2 


<3 


HD 155802 


s 


34.27 


5010±15 


4.46±0.04 


-0.15±0.02 


13±4 


0.96±0.13 


2 


5±1 


HD 171067 


s 


39.73 


5520±35 


4.03±0.04 


-0.15±0.03 


12±4 


1.43±0.16 


3 


<3 


HD 184591 


s 


4.58 


4800±40 


2.82±0.07 


-0.18±0.03 


14±6 


0.80±0.24 


3 


4.7±0.8 


HD 218739 


s 


34.06 


5670±20 


4.19±0.05 


-0.07±0.02 


135±6 


2.80±0.05 


3 


7.6±0.2 


HD 40891 


a 


32.95 


5140±25 


4.37±0.05 


-0.10±0.02 


11±4 


1.03±0.16 


3 




HD 62668 


a 


4.97 


4660±40 


2.92±0.10 


-0.30±0.04 


82±8 


1.60±0.10 


3 


21±1 


HD 66553 


a 


28.03 


5275±40 


4.33±0.04 


+0.09±0.02 


12±7 


1.20±0.27 


2 




HD 82159 


a 


21.11 


5150±65 


4.62±0.08 


-0.13±0.03 


109±7 


2.22±0.10 


2 


12.5±0.5 


HD 82841 


a 


4.09 


4620±30 


2.94±0.07 


-0.30±0.03 


19±6 


0.74±0.19 


3 


8.1±0.3 


HIP 50072 


a 


3.15 


4615±39 


3.08±0.08 


-0.33±0.04 


22±6 


0.80±0.21 


3 


20.8±0.8 


HD 112099 


a 


38.12 


5095±20 


4.37±0.04 


-0.14±0.02 


11±4 


1.00±0.16 


2 


4.1±0.2 


HIP 63322 


a 


26.24 


4995±27 


4.48±0.07 


-0.23±0.04 


142±8 


2.19±0.07 


2 


6±1 


HIP 63442 


a 


3.31 


4530±30 


2.85±0.08 


-0.35±0.04 


41±6 


1.07±0.13 


3 


22.3±0.6 


HD 138157 


a 


5.07 


4845±82 


3.24±0.12 


-0.36±0.06 


20±5 


1.03±0.22 


3 


40±1 


HD 147866 


a 


5.18 


4580±15 


2.71±0.05 


-0.24±0.02 


13±4 


0.49±0.14 


3 


5±1 


HD 150202 


a 


3.75 


5010±25 


3.06±0.07 


-0.12±0.03 


25±5 


1.32±0.16 


3 


7.8±0.2 


eUMi 


a 


9.41 


5215±50 


3.21±0.08 


-0.25±0.04 


38±6 


1.80±0.16 


3 


25.6±0.2 


HD 190642 


a 


4.48 


4760±65 


2.96±0.11 


-0.47±0.05 


19±6 


0.90±0.23 


3 


19.6±1.0 


HD 202109 


a 


21.62 


4910±15 


2.14±0.04 


-0.08±0.02 


21±5 


1.13±0.15 


3 




HD553 


a 


4.74 


4645±90 


3.38±0.36 


-0.23±0.08 


22±4 


0.85±0.22 


3 


41±2 








5940+1 50 


4.40±0.77 


-0 4 1 ±0 26 


21±6 


2 n7±0 26 


2 


27± 


HIP 999 


a 


24.69 


5280±70 


4.40±0.30 


-0.31 ±0.04 


20±4 


1.44±0.19 


2 


27±1 




b 




4840±150 


3.6±1.2 


+0.3±0.3 


15±8 


0.83±0.42 


3 


50±2 


HD 8997 


a 


43.16 


5060±65 


4.53±0.25 


-0.12±0.04 


10±6 


0.91±0.27 


2 


3.0±0.5 




b 




4525±50 


4.74±0.26 


+0.02±0.18 


11±4 


0.34±0.21 


3 


8±1 


HD 9902 


a 




5130±80 


4.0±0.4 


-0.50±0.16 


52±7 


1.85±0.16 


2 


12±1 




b 




5715±150 


3.37±0.20 


-0.38±0.15 


26±8 


1.97±0.31 


5 


2.0±0.5 


HD 16884 


a 




4500±50 


2.4±0.4 


+0.29±0.11 


31±7 


0.85±0.20 


3 


6±1 




c 




4535±70 


3.4±0.5 


-0.26±0.15 


10±9 


0.36±0.38 


3 


44±2 


HD 18955 


a 


20.56 


5410±80 


4.84±0.10 


+0.03±0.06 


9±2 


1.25±0.13 


2 


9±1 




b 




5095±150 


4.8±0.3 


+0.22±0.17 


8±2 


0.90±0.23 


2 


5±1 


SAO 15 1224 


a 




4595±110 


3.7±0.4 


-0.18±0.14 


45±7 


1.20±0.23 


3 


45±2 




b 




5420±105 


4.8±0.3 


-0.09±0.12 


22±6 


1.67±0.25 


2 


6±2 


HD 45762 


a 


3.45 


4630±100 


3.0±0.7 


-0.48±0.24 


31±5 


1.02±0.22 


3 


60±2 




c 




6015±150 


3.5±1.0 


-0.19±0.18 






5 


16±2 


HD 50255 


a 


30.26 


5580±100 


4.07±0.12 


-0.18±0.11 


45±6 


2.19±0.18 


2 


3.0±1.0 




b 




4980±150 


3.7±1.0 


-1.0±0.8 


9±4 


0.81±0.30 


3 


5±1 
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Table 7 (continued) 



Name 


CD 


TT 

(mas) 


J- eft 

(K) 


logs 
(ms"^) 


[te/H\ 


(mA) 


log n[Li) 




V sin i 
(km s^^) 


HD 61994 


a 


35.13 


5630±150 


4.13±0.11 


+0.03±0.12 


17±4 


1.69±0.23 


2 


4±1 




b 




4775±150 


4.6±0.6 


^0.27±0.33 


9±4 


0.56±0.33 


3 


9±1 


HD 73512 


a 


39.28 


5110±115 


4.58±0.12 


-0.25±0.04 


11±5 


1.00±0.17 


2 


5±1 




b 




4600±130 


4.4±0.5 


-0.37±0.17 


11±3 


0.41±0.27 


3 


5±1 


HD 82286 


a 


9.57 


4800±40 


3.73±0.31 


-0.05±0.09 


21±5 


1.00±0.17 


3 


37.5±1.5 




b 




4785±150 


3.92±0.35 


-0.42±0.20 


24±6 


1.04±0.32 


3 


40.5±2 


HD 93915 


a 


22.78 


5520±105 


4.28±0.24 


-0.29±0.09 


12±4 


1.43±0.22 


2 


5±1 




b 




5315±120 


4.25±0.24 


-0.22±0.09 


9±3 


1.17±0.22 


2 


3±1 


HD 95559 


a 


18.43 


5090±110 


4.47±0.27 


-0.10±0.12 


16±6 


1.14±0.28 


2 


33±2 




b 




5065±150 


4.56±0.33 


+0.07±0.12 


20±6 


1.23±0.32 


2 


33±2 


HD 105575 


a 


19.77 


5650±150 


4.73±0.11 


-0.46±0.08 


14±2 


1.61±0.17 


2 


42±15 


HD 106855 


a 


22.88 


4750±110 


4.48±0.33 


-0.87±0.19 


22±4 


0.91±0.24 


2 


30±2 




b 




4950±150 


4.0±0.8 


-0.5±0.4 


11±4 


0.85±0.30 


3 


5±2 


HD 109011 


a 


Ai.n 


5030±75 


4.59±0.27 


-0.17±0.10 


30±5 


1.42±0.19 


2 


5±1 




b 




4900±150 


3.7±1.0 


+0.2±0.3 


13±5 


0.84±0.33 


3 


6±1 


HD 111487 


a 




5500±150 


4.47±0.28 


-0.36±0.11 


13±4 


1.45±0.26 


2 


42±2 




b 




4700±150 


4.6±0.5 


-0.1 ±0.5 


14±6 


0.64±0.38 


3 


28±2 


HD 112859 


a 


5.24 


4780±120 


3.41±0.38 


-0.24±0.15 


47±7 


1.45±0.23 


3 


20±1 




b 




6150±150 


2.5±0.5 


+0.05±0.3 






5 


5±1 


HD 127068 


a 


9.75 


4850±135 


3.61±0.27 


-0.48±0.07 


76±5 


1.79±0.19 


3 


6±1 




b 




5560±150 


3.9±0.8 


-0.39±0.14 


10±5 


1.46±0.31 


3 


8±2 


HIP 77210 


a 


20.73 


5190±120 


4.41±0.26 


-0.43±0.10 


14±4 


1.18±0.25 


2 


3±1 




b 




4720±150 


4.74±0.31 


+0.08±0.11 


16±5 


0.73±0.33 


3 


5±1 


HD 142680 


a 


28.23 


4885±110 


4.35±0.26 


-0.22±0.05 


5±3 


0.54±0.24 


3 


3±1 




b 




4500±150 


4.6±0.6 


-0.51±0.36 


8±3 


0.21±0.29 


3 


5±1 


HD 143937 


a 


23.73 


5160±140 


4.48±0.16 


-0.30±0.13 


12±6 


1.11±0.33 


2 


55±2 




b 




4730±150 


4.72±0.17 


-0.01±0.13 


13±6 


0.65±0.35 


3 


45±2 


HD 197913 


a 


17±6° 


5520±100 


4.18±0.32 


-0.06±0.11 


35±4 


2.00±0.17 


2 


5±1 




b 




5340±120 


4.20±0.45 


-0.08±0.10 


19±6 


1.47±0.28 


2 


6±1 


HD 199967 


a 


13.57 


5850±100 


4.26±0.22 


-0.24±0.17 


28±7 


2.13±0.22 


2 


11±1 




b 




5835±150 


4.4±0.4 


-0.25±0.13 


19±8 


1.91±0.33 


2 


9±1 


HD 226099 


a 


35.46 


5200±80 


4.45±0.34 


-0.92±0.23 


12±4 


1.12±0.21 


2 


2±1 




b 




4560±80 


4.20±0.36 


-l.liO.l 


9±4 


0.30±0.25 


3 


4±1 


HD 237944 


a 
b 
c 


10.62 


(5630)'' 
(5560)'' 
(4900)'' 










3 
3 
3 


11±2 

9±2 

11±3 



"from Jenkins {1952), 

''assumed from spectral classification. 



transformation and the effective temperatures and gravities 
from Table |7] 



4.6 Ha- core flux 

We measure the central 1-A portion of the Ha line from 
our continuum-normalized spectra and relate it to the abso- 
lute continuum flux, Jc, expected at Ha. The latter is ob- 
tained from the relations provided by Hall ( 1996| l for vari- 
ous Morgan-Keenan (MK) classes and color ranges, i.e. 



log J"c = 7.538 - 1.081 {B - V) 
for MK I - V and < (B - 1/) < 1.4 , 
log J"c = 7.518 - 1.236 {V - R) 
for MK V and < (y - i?) < 1.4 , 



(1) 



(2) 



logJ-c^ 7.576- 1.447 (y-i?) (3) 
for MK I - IV and < (T^ - i?) < 1.8 . 

The stellar surface Ha-core flux in ergcm~^s~^ is then 
computed from the measured 1-A equivalent width, Wcore, 
under the spectrum and zero intensity. 



(4) 



In this way the equivalent width increases with increased 
line-core filling due to chromospheric emission while the 
continuum flux reflects the dependency on stellar radius 
and Tcff ■ Rotational line broadening is accounted for by ap- 
plying a similar correction to W^corc as done by Pasquini 
& Pallavicini ( |1991| l. We obtain these corrections by artifi- 
cially broadening the spectrum of HD 153525, a KO dwarf 
with an otherwise unresolved usini, in steps of 5 kms~^ 
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Fig. 5 A comparison of absolute continuum flux at Ha from different calibrations. Only our single and SBl stars are 
plotted, a Fluxes based on the effective temperatures from Table |7] converted to B — V with the T^s vs. B ~ V calibration 
of Flower ( 1996| ) versus the flux based on the B — V listed in the Tycho-2 catalog, b Again our effective temperatures 
converted to B — V but with the Toff calibration of Gray (2005 ) versus the B — V listed in the Tycho-2 catalog, c Our 
effective temperatures converted to y — i? with the Tcff calibration of Gray ( 2005 | l versus the B — V"-based fluxes listed in 
the Tycho-2 catalog. Minimized rms errors of the flux residuals are achieved in panel a. The dashed line is just a 45-degree 
line. The solid lines are non-weighted regression fits. 



up to 60 km s^^ and then measure the Ha line as done pre- 
viously. The correction factor, A, for each stellar measure- 
ment is obtained from the fit to these simulations by insert- 
ing the stellar values for v sin i from Table |7| 

A = A + B vsini + C vsinf + D vsini^ . (5) 

The coefficients are yl=-0.00068990582, 
5=0.00029478396, (7=0.00011360807, 7:>=-9.4974126e- 
07. Note that the corrections are negligible for 
vsini <10 kms^^ but reach a maximum for HD 45762 
(vsini=60 kms-i) of 30%. Table [s] lists the corrected 
fluxes. 

The line-core fluxes from Eq. |4] are not corrected for 
photospheric light. Such a correction for the Ha line core 
is not without problems because several line-formation pro- 
cesses are at work at the same time and significant star-to- 



with the color transformation from Flower ( |1996| l. Two sec- 
ondary choices were fluxes from our T^g values but con- 
verted to B -V eiudV - R using the tables in Gray (|2005 i. 



star differences are expected (see Cram & Mullan 1985 1. 



There is no readily available photospheric Ha-core flux tab- 
ulated for standard model atmospheres. Therefore, previous 
attempts relied on empirical methods, e.g., the subtraction 
of an observed profile of a bona-fide inactive star of other- 



wise similar quantities (e.g. Soderblom et al. 1993b I or the 
subtraction of an empirically predetermined minimum flux 



from a large sample of stars (e.g. Cincunegui et al. 2007 1. 
However, as was shown nicely in the work by Cincunegui et 
al. ( |2007[ ), this photospheric correction significantly wors- 
ened the relationship with Ca II H&K flux. The conclusion 
is that empirical photospheric corrections may not be fully 
photospheric in origin. In the present paper, we refrained 
from any photospheric correction and give combined pho- 
tospheric H-chromospheric fluxes as observed. 

The largest contribution to the flux error is due to the 
continuum-flux calibration and comes from systematic dif- 
ferences in the Tcff -color relations and the error in Toff- We 
compared fluxes based on four different relations. Our two 
prime choices for are based, firstly, on the B — V from 
the Tycho-2 catalog (H0g et al. |2000| l and, secondly, on the 
B — V from the conversion of our measured T^ff in Table |7] 



The inter comparison in Fig. |5] shows minimized rms er- 
rors of the flux residuals for the two prime choices, i.e. the 
Tycho-2 B — V^-based fluxes vs. the fluxes from measured 
Tcff 's with the color conversion from Flower. It gave an rms 
of 0.45 in units of 10^ ergcm^^s^^. Comparable rms is 
obtained for the Tycho-2 B — V^-based fluxes vs. Gray's 
B — V based fluxes but significantly larger (0.75) if Gray's 
V — R based fluxes are used instead. The calibration rms 
errors per star are listed in Table [8] along with the line-core 
fluxes based on the Tycho-2 B — V color for the single and 
SBl stars, and on flie Flower ^199^ T^g-{B - V) trans- 
formation for the SB2s using the conversion in Eq. [T] These 
rms values are useful when judging the reality of time varia- 
tions, if present. An exceptionally large rms value indicates 
most likely that the calibration assumptions in Eq.[T]-3 were 
not fully appropriate, e.g. for SAO150676, which is likely a 
pre-main sequence star where the applicability of the T^g- 
(B — V) conversion is questionable. 

The number of spectra used for the Ha analysis is 
smaller than for the radial-velocity curves because we use 
only the ones with the best S/N and, in case of the SB2 
binaries, when the line separation is large enough. We mea- 
sure single-lined spectra just like single stars, which have 
continua presumably unaffected by the companion. The few 
systems where we see a Ha line from the secondary com- 
ponent, but not its photospheric lines, were also treated like 
single stars. The regular SB2 spectra, on the other hand, 
were first multiplied by the respective photospheric light ra- 
tio, i?(Fel). It is obtained from line pairs near Ha (with 
highest weight given to the unblended Fe I 654.6 nm line) 
before measuring the Ha flux for each component. These 
fluxes are generally less precise because the available S/N 
ratio in the continuum is distributed among the two com- 
ponents. Table |8] lists the measurements for all single stars 



© 20 1 2 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN (2012) 



23 



and single-lined binary stars as well as for all individual 
stars in the double-lined binaries. The triple-lined system 
HD 105575 was too uncertain to be measured reliably and 
only the primary is included in Table [8] 

5 Notes to individual systems 

HD 553 = V741 Cas. The system was discovered to be an 
eclipsing binary by the Hipparcos satellite. In paper I, we 
had verified its eclipsing nature and shown that the system 
is double lined. A first orbit was published shortly thereafter 
by Duemmler et al. ( 2002| followed by an independent or- 
bit by Griffin ( 2003| l. Both their orbits were circular with 
a period of 9.0599 d and a solution rms residual of 0.56 
and 1.36 kms^^, respectively. The component's rms in the 
Duemmler et al. orbit were 0.59 and 0.97 km s"^ Our new 
orbit is also assumed to be circular The forced e=0 solution 
gives residuals of 0.41 and 0.39 kms~^ for the two com- 
ponents, respectively. There is no conclusive period in the 
Ha-core flux but a peak at 4.5 d could be interpreted to be 
~ Porb/2. The out-of-eclipse light variations have a V am- 
plitude of 0.10 mag and suggest a photometric period of 
9.061 d (eclipse depth is 0.28 mag in V). 

HIP 999 = LN Peg. An initial SBl orbit was determined 
by Latham et al. ( |1988| l while Fekel et al. ( |1999[ ) presented 
a SB2 orbit and found the star to be at least a triple system. 
Our radial-velocity residuals indicate a long-term period of 
around 4.2 years with a full amplitude of 12 kms^^zero 
eccentricity (see Fig. l2|. This verifies the preference put 
forward by Fekel et al p999| l that the long-period orbit is 
rather 4 years with zero eccentricity than 8 years with 0.5 
eccentricity. The long orbit of 1555 d is listed in Table|4]and 
was obtained by including the Fekel et al. data with a shift of 
H-0.50 kms~^. The total time span of data was then 9430 d 
with 213 velocities. The STELLA data alone would give el- 
ements of ii:3=5.654±0.085 kms-\ To=2,454,117.4±8.5, 
P=1310±16 d. The Ha line from the secondary star in the 
close pair appears in emission. The third component c is 
not directly seen in the spectra but is recognized only due 
to the systematic radial velocity zero-point variations of the 
ab pair The a6-pair orbit is given in Table |5] and shown in 
Fig. [3] 

HD 8997 = EO Psc. Griffin's ( fT987] l SB2-orbit is basi- 
cally confirmed but now considerably more precise. Our low 
rms of wO.4 km s^^ for the a and b component radial veloc- 
ities allowed to search for residual variations and a period of 
88.3±0.5 d is found for the primary with a full amplitude of 
0.20 km s~^. No period is seen in the secondary-star residu- 
als. If above period would be interpreted to be the rotational 
period of the primary, and combined with our v sin i mea- 
sures of 3 kms^^, the minimum radius would be 5.2 Rq. 
This value is clearly in disagreement with the distance of 
23 pc and a system brightness of 7.74 mag and B — V^=0.96. 
This period must be of other origin or is even spurious. The 
photometry yields a period of 10.54±0.11 d, with a small 
amplitude of 0.02 mag in y. This period suggests i?sini 



of «0.6 R0, in good agreement with the distance and the 
spectrum synthesis that indicates Tcfr=5060±65 K, a log g 
of 4.5, and a metallicity of -0. 12±0.04, suggesting a K2 
dwarf star. 

HD 9902 = BG Psc. This system contains two almost 
equally massive stars with a mass ratio of 0.9750±0.0026 
in a 25.3722-d orbit with e=0.5068. We refer to the slightly 
more massive and more rapidly rotating star as the primary. 
Note that the rms for the orbit for the secondary (154ms~^) 
is lower than for the primary due to its sharper lines. Our 
photometry of the combined system reveals a period of 
7.44±0.02 d that we interpret to be the rotation period of 
the primary (usini=12 kms^^) rather than that of the sec- 
ondary (vsmi=2 kms^^). Also note that the V amplitude 
of 0'.^05 in Table|6]is for the combined light and not for, say, 
the primary. The expected pseudo-synchronized rotation pe- 
riods are w8.9 d according to Hut ( 1981| l. The primary is 



thus likely rotating faster than synchronous. 

HD 16884. The star is double lined from the cross- 
correlation functions but its spectral lines do not ap- 
pear to originate from a single physical pair. Most likely 
HD 16884 is a quadruple system and we see lines from 
the two primaries a and c, with each system being sin- 
gle lined, i.e., SBl-nSBl. The orbital periods for the two 
pairs (each 106.63 d) appear to be identical within its er- 
rors, and so does the systemic velocity (h-2.2 kms~^). The 
K-amplitude of the a and c components are 30.0 km s^^ and 
16.3 km s^^, respectively, while their v sin i's are 6 km s~^ 
and «44 kms^^, respectively, which means that the radial 
velocity attributed to the c component can not be caused 
by rotationally modulated line-profile variations of the a 
component, its v sin i is too small. Not much can be said 
about the orbital period ab-cd but the two nodal lines ap- 
pear co-aligned and did not change recognizably during our 
observations, its orbital period must be very long. Griffin 
2009| l obtained a sparsely sampled single-lined orbit from 



14 of his own velocities, an orbit which is superseded in the 
present paper. We also note that the stronger-lined a star's 
Ha line profile appears as a broad and asymmetric emission 
line (full width at continuum of up to 400 kms^^) with a 
deep central absorption that follows the radial velocities of 
the a component (Fig. |6]). Our photometry reveals a period 
of 65.44±0.05 d and an amplitude of 0.16 mag in V for 
the combined stars. It is likely that this is the rotation pe- 
riod of the brighter a component, a K4 giant according to 
Griffin ^ 2009) 1. His classification is in rough agreement with 
our Tcff of 4500 K and log g of 2.4±0.4 from the spectrum 
synthesis. With a iisini of 6 kms~^ its minimum radius, 
i?sini, would be 7.8±1.3 R©. However, a K4 giant is ex- 
pected to have a nominal radius of approximately 25 R©. 
Thus, if indeed a K4 giant, i would be only « 18° and to- 
gether with an assumed mass of 2.3 for the primary, 
this contradicts the large mass function of 0.2972±0.0013, 
which would otherwise imply an unreasonable secondary 
mass of order 7.5 M©. Doubling the photometric period 
to get the true rotation period would lead to an inclination 
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HD 16884 




0.4^ 



6555 ' 6560 6565 6570 

Fig. 6 Ha profile variations of HD 16884 in 2007/08. Spectra are corrected for orbital radial velocity of the primary. The 
photospheric absorption lines belong to the primary of the ab system, a K giant, and agree with the radial velocities from 
Ha. It is save to conclude that the majority of the Ha profile originates from this K giant. The signature from the primary 
component of the cd system, i.e. the secondary in this composite spectrum, is lost in this representation. See text. 



much closer to 90° and, lowering the assumed primary mass 
to 2.0 Mq the secondary mass gets closer to a minimum 
mass of ~ 1.6 Mq, still rather high for an "unseen" star. Cer- 
tainly, this system warrants further investigation. 

HD 18645 = FU Cet. Single star. Our photometry sug- 
gests a period of 21 .5 d that we interpret to be its rotation pe- 
riod. A residual period of 3.069 d is seen from the radial ve- 
locities but with a very weak amplitude of 0. 1 km s^ ^ . Addi- 
tionally, the Ha-core emission shows a seasonal trend and, 
after its removal, a period of «17 d. Although this period is 
badly determined and inconclusive it may suggest that the 
photometric period of 21 .5 d is to be preferred against the 3- 
d period. Our vsini measure of 12.0±0.2 kms^^ converts 
then to a minimum stellar radius of 5.1 R©, appropriate for 
its mid G giant classification. 

HD 18955 = IR Eri. This system consists of two 
dwarf stars of mass ratio 0.88 in a 43-day orbit with a 
high eccentricity of 0.7587±0.0005. The expected pseudo- 
synchronization period would be a mere 1 1% of the orbital 
period, i.e. ?s5 d. Unfortunately, our photometry did not al- 
low to obtain a period despite that the data indicate variabil- 
ity at the 2.3(7 level. Both components show a weak resid- 
ual radial-velocity variation with a period of 25-26 d and 
an amplitude of 0.5 kms^^. This could be the stellar ro- 
tation period but the data are too noisy to be conclusive. 
It appears that the secondary is cooler than the primary by 
w300 K, and also slower rotating (wsini=5 kms^^ vs. 9, 
respectively). 

HD 23551 = MM Cam. The star appears to be single. 
A radial velocity of -7.5±0.3 kms~^ was reported in the 
Pulkova survey by Gontcharov ( |2006| in good agreement 
with the average -7.38 ± 0.22 by Famaey et al. ( |2005| ). Our 
two velocities from KPNO in paper I were —9.5 from the 
blue spectrum and — 7.4 ± 0.2 from the red spectrum, also 
in reasonable agreement. A peak with a period of «80 d is 
existent in the Ha flux and the photometry suggests a clear 
period of 8 1 .69±0.03 d that we interpret to be the rotation 
period of the star. Note that we observed the star with both 



APTs (T6 and T7) and that the independent photometric pe- 
riod from T7 data is 81.2±0.5 d, in good agreement with 
the more precise T6 period. Its vsini of 6.4±0.2 kms^^ 
converts then to a minimum radius of 10.3 Rq in agreement 
with log 5=3.0 from our spectrum synthesis. We also note its 
significant metal deficiency of -0. 1 1±0.02. No Li I 6707.8- 
A absorption feature is seen above our detection threshold 
of wlOmA (Fig.^. 

HD 24053. This star first appeared in the literature in 
the uvby calibration for FGK stars by Olsen, e.g., Olsen 
I994| l. Its effective temperature was determined spectro- 



scopically from spectral line ratios by Kovtyukh et al. 
( 2003| l to 5723 K. The revised Geneva-Copenhagen cat- 
alog (Holmberg et al. '2007) lists it with a metallicity of 
—0.04 and 5640 K and suggests an age of 0.2 Gyr. From 
2-MASS colors Masana et al. ( 12006) 1 derived 5711 K and 
the N2K-consortium (Robinson et al. |2007| ) derived 5699 K, 
a metallicity of +0.09, and log (7=4.45 cms~^. Earlier ve- 
locities were given in our paper I (5.3±1.4 kms^^ from 
a blue spectrum and 2.9±0.6 kms^^ from a red spectrum 
two days apart) while the catalog by Gontcharov (| 2006[ l 
lists 4.1 ±0.2 kms^^. Two spectra are given in the Elodie 
archive. We have 77 STELLA spectra and all of them ap- 
pear to be constant. The average values from our spectrum 
synthesis yields 5600±25 K, log5=4.09±0.03, and a metal- 
licity of -0.01±0.02. No significant variations are seen in 
Ha. 



SAO 150676 = AILep. Cutispoto et al. ( fT999) ) fisted the 
star in their ROSAT follow-up observations. They noted 
a spectral classification of G2V but added "(PMS?)" in- 
dicating possibly pre main sequence. In the same study 
a photometric variability of 0™02 in V with a period 
of 1.78 d was reported. The authors also measured a 
vsini of 28±2 kms^^ in reasonable agreement with 
25.6±2 km s^^ from our paper I . Their two radial velocity 
measures agreed to within their errors and were averaged 
to H-25.8±1.6 kms^^. In paper I we had given velocities 
of H-26.1±1.6 kms^^ from a single red- wavelength spec- 
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Fig. 7 Li I 6707.8- A line region for the KOIII component 
of the SB-1 binary HD 62668 (top spectrum; shifted by 
+0.2) and for the single KOIII star HD 23551 (bottom spec- 
trum). Notice the broad Li I absorption line in HD 62668 
and its absence in HD 23551. 



trumandH-21.5±5.3 kms^^ and — 18.5±2.1 kms^^ for two 
components from a single Ca II H& K emission-line spec- 
trum. The latter must have been an erroneous measurement 
because it is not verified in the present paper In paper I, we 
also measured the large equivalent width of 195mA from the 
Li I 670.8-nm line, in agreement with the suggestion of be- 
ing a pre-main sequence star Our 93 STELLA radial veloci- 
ties have an unreasonable high scatter of 0.53 km s^^ which 
is likely due to rotational modulation of a spotted surface. 
A Lomb periodogram shows two equal peaks at periods of 
0.64 d and L79 d. Unfortunately, the Ha fluxes do not show 
the same peaks but favor a (poorly determined) 1 .2-day pe- 
riod. However, our photometry shows a clear 1 .79-d period, 
in good agreement with Cutispoto et al. (1999) and one of 
the two radial-velocity periods. We interpret it to be the ro- 
tation period of the star. 

SAO 151224 = HY CMa. The star is an echpsing SB2 
binary in a circular 4.97769 1-d orbit. Photometry from late 



1993 (Cutispoto et al. |1995| [T999) revealed the echpsing 
nature while follow-up photometry in 1994 by Cutispoto 



et al.(2003l failed to confirm the eclipses and the author's 
stated to the contrary that "No primary eclipses are appar- 
ent." Our photometry verifies the eclipsing nature of the 
light variations with a period of 4.98 d and an eclipse depth 
of 0.68 mag in V . The residual radial velocities of the pri- 
mary star indicate a period of 5.25 d, likely due to rota- 
tion of its spotted surface which partly explains the unex- 
pected high rms of its orbital fit of 1.13 kms^^. The vsmi 
of the primary of 45 kms^^ converts directly into a radius 
of 4.4 Rq. 

HD 40891. We see no traces of Call H&K emission in 
our echelle spectra. Therefore, its emission can be classi- 
fied at best "very weak". Note that in paper I, we listed the 



star in the table for stars with emission, but did not show 
a plot of the spectrum. This might indicate a logistics error 
in our paper I. The star is an SBl in an eccentric orbit with 
e=0.5105±0.0007 and a period of 49.0823 d. The expected 
pseudo-synchronization period would be wl7 d. 

HD 45762. Griffin pOTOl l noted in his Table 1 "Bi- 
nary, P^60 days" but otherwise gave no information. We 
find the star to be a double-lined triple system with the 
more massive component being an SBl with a period of 
wl/3 that of the wider c component. The (ab) — c orbit 
has a period of 59.9363±0.0016 d and a small but non- 
zero eccentricity of 0.0150±0.0009. Its rms of the residu- 
als is 0.533 kms^^ which increases to 0.88 kms^^ if an 
e=0 solution is enforced. The mass ratio c/(a + 6) is 0.23 
from an averaged ac solution (shown as a dotted line in 
the respective panel in Fig. [3]l. The close SBl orbit has 
only a small amplitude and we find the following pre- 
liminary elements (e=0 assumed); K=l .1%±.Q.11 kms^^, 
ro=2,454,116.84±0.49, and P(a6)=19.933±0.014 d. Note 
that this star's vsmi is 60 kms^^, as compared to the c- 
component's 16 kms^^, and is likely an asynchronous ro- 
tator because a Piot=19.933 d would suggest a large mini- 
mum radius of 23.6 R©. Our photometry shows a double- 
humped light curve with a clear 62-day periodicity but must 
be partly due to an ellipticity effect because of its phase co- 
herence with the orbit. The c component has an expected 
minimum radius of 19.6 R© if Prot=62 d. Obviously, both 
stars are cool giants, as already indicated by the small Hip- 
parcos parallax and the y=8.3-mag brightness. 

HD 50255. The star is a SB2 in an eccentric 455-d orbit. 
The scatter of the residual radial velocities of the secondary 
component (0.167 kms^^) is larger than expected for a 
star with vsuii=5 kms~^. A Lomb periodogram shows 
significant excess power at 485 d with an amplitude of 
0.65 kms~^, close to the orbital period of 455 d, but given 
the comparable short coverage of 1471 d of the entire data 
set it could be also a sampling effect and remains to be con- 
firmed. Tentatively, we interpret it to be the rotation period 
of the secondary. Note that our photometry did not reveal 
any clear Ught variations above, say, 0.05 mag. 

HD 61994. There are 102 secondary-star velocities de- 
tected from the STELLA cross-correlation functions and an 
eccentric SB2 orbit is given (P=552.8 d, e=0.423). Duquen- 
noy & Mayor (1988) had discovered its binarity and also 
presented a (marginal) SB2 orbit based on 9 CORAVEL 
measurements. The visual brightness difference between the 
two components of 2.6 mag is close to its detection limit. 
The average line ratio in our spectra is 8.5. The expected 
pseudo-synchronized rotation period is w250 d, but neither 
the velocity residuals nor the Ha-core flux show traces of 
any significant period. 

HD 62668 = BM Lyn. Our orbit is assumed circu- 
lar with a period of 69.3023 d. A formal eccentricity of 
0.0043±0.001 1 is judged too uncertain. The photometry 
is best fit with a period of 67.470±0.007 d, close but not 
identical to the orbital period. The Ha absorption line also 



www.an-joumal.org 



© 2012 WILEY- VCH Verlag GmbH & Co. KGaA, Weinhdm 



26 



K. G. Strassmeier, M. Weber, T. Granzer, S. Jarvinen: Rotation, activity, and lithium abundance in cool binary stars 



Table 8 Ha absolute emission-line fluxes. Column SB denotes the binary component in case the system is SB2 or 
even SB3 (S means single star, a means primary component, b secondary component, c tertiary component). (B-V) and 
(V-R)c are the expected Johnson-Cousins indices from our Toff and the transformation from Flower ( |1996| l and Gray 
( 2005 1, respectively, that we adopted for the flux calibration. N is the number of spectra used. Wcorc is the equivalent 
width in mA in the central 1-A bandpass and its internal measuring error. (J^(Ha)) is the average 1-A Ha-core flux in 
10^ erg cm^^s^^ and its propagated error from the equivalent width measure plus the flux calibration, "rms" is its standard 
deviation from the four different flux calibrations for comparison. ST is the rms residual of J^(Ha) over the time span of 
our entire observations. i?(Fe l) is the photospheric line ratio secondary/primary from Fe I 654.6 nm that is used to weight 
the components' Ha equivalent widths for SB2 systems. 



Name 


SB 


N 


B-V 


V-R 
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(-F(Ha)) 
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(mag) 


(mA) 


(xlO^ 


erg cm" 




(-) 


HD 18645 


S 


64 
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Table 8 (continued) 



Name SB TV B-V V-R Wcore (J^iUa)) rms ST(Ha) i?(Fel) 

(mag) (mA) (x 10^ ergcm^^s"-^) (-) 
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shows a clear time-dependent variability (Fig. |8]l but no 
period. The residual radial velocities indicate a period of 
22.35 d, which we believe is an alias of the 69-d period. A 
weak lithium line at 6707.8 A is seen (a representative plot 
is shown in Fig.]?}. 

HD 66553. There are some secondary-star velocities de- 
tected from our cross-correlation functions but close-up in- 
spection of the spectra did not unambiguously show spec- 
tral lines from it. Therefore, we refrained from giving an 
SB2 orbit. This system shows the most eccentric orbit in 
our sample (e=0.8719±0.0028) at a period of 75.8996 d. 



Latham et al. (20021 discovered the radial-velocity varia- 



tions and determined an SBl orbit. "The ninth catalogue of 
spectroscopic binary orbits" (Pourbaix et al. 2004) lists an- 
other independent SB-1 orbit from Imbert ( (2006 j ), in agree- 
ment with the Latham et al. orbit. The expected pseudo- 
synchronization period would be ?s3.4 d, i.e. just «4% of 



the orbital period. No signs of periodic velocity residuals or 
Ha-flux variations are evident. Thus, the rotation period of 
the star still remains to be determined. Our spectrum synthe- 
sis gives roff=5275±41 K, log5=4.33±0.04, and a higher 
than usual metallicity of H-0.09±0.02. 

HD 73512. A SB2 orbit with a period of 128.241 d and 
e=0.26240±0.00045 is given. The mass ratio is 0.88 while 
the component's line-depth ratio is 2.2 with both stars hav- 
ing the same vsini of «5 kms^^. The expected pseudo- 
synchronized rotation period is w90 d. The star appears con- 
stant in our photometry. 

HD 76799. The star shows radial-velocity variations 
with an amplitude of 0.23 kms^^, peak to valley, with 
an uncertainty of ±0.04 kms^^ if phased with the most- 
likely period of 50.13±0.15 d from a Lomb-Scargle analy- 
sis. A periodogram from four years of Ha line fluxes shows 
two peaks with periods of «64 d and 55 d. Although not 
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Fig.8 Ha profile variations of, left, HD 62668 (SBl) in 2007 and, right, HD 82286 (SB2) in 2010. Spectra are corrected 
for orbital radial velocity. Note that the spectral lines from the secondary star in HD 82286 appear as apparently increased 
scatter around the primary-star lines while the sharp lines are actually water-vapor lines. 



quite in agreement with the photospheric radial-velocity pe- 
riod, they are in the same range. We interpret both varia- 
tions to be due to stellar rotation. The star is otherwise a 
single star. Our spectrum synthesis gave rcff =4915±24 K, 
log (/=3.26±0.07, and a significantly sub-solar metallicity 
of -0.12±0.03. The usini of 7.8±0.3 kms^^ suggests a 
minimum radius of 7.7 Rq, in agreement with the low grav- 
ity and the small parallax. 

HIP 46634 = BD+11 2052 B. This star is the B compo- 
nent of the visual binary ADS 7406, about 14" away from 
the brighter A component. Galvez et al. ( |2005 i found it to 
exhibit a constant radial velocity of H-27.07±0.04 kms^^ 
within their time of observation in 2004. The A component 
is HD 82159 and was identified in our KPNO H&K sur- 
vey as an active single-lined spectroscopic binary. Our cur- 
rent STELLA/SES spectra show the B component to have 
a constant radial velocity of H-27.927±0.047 kms~^. It is 
apparently single. There are two equally strong peaks of 
372 d and 185 d in the radial velocities with an amplitude 
of 0.2 kms^^. The Ha-line fluxes show peaks at 120 d and 
298 d indicating some long-term variations. None of the pe- 
riods is conclusive though and the true rotation period re- 
mains to be determined. 

HD 82159 = GSLeo = BD+11 2052 A. This is the bright 
component of the double star ADS 7406, denoted compo- 
nent A. Its orbital elements with a period of 4.8122 days 
and e=0.176 were given by Galvez et al. ( |2005 20061 from 
data taken within nine nights in March/ April 2004. Griffin 
& Filiz Ak (2010) had revised this orbit from 30 new Cam- 
bridge CORAVEL measurements from 2008/09. Our SBl 
orbit is computed from 107 STELLA measures spread over 
four years with a period of 3.855871±0.0000066 d and an 



eccentricity of 0.2640±0.0012 and supersedes yet the Grif- 
fin orbit. The velocity residuals are systematic and indicate 
a period of 3.50 d, significantly longer than our photometric 
period of 3.054 d. Note that our photometry is polluted by 
the B component HIP 46634 which, however, should not af- 
fect the period determination. The Ha-core emission fluxes 
indicate, besides an annual trend that must be removed, a 
residual period of 2.97 d, in good agreement with the pho- 
tometric period. Therefore, we feel confident to interpret the 
photometric period to be the rotation period of the primary 
star. Its expected pseudo-synchronization period is «2.7 d. 
Obviously, the system did not reach equilibrium yet and the 
primary is still rotating asynchronously. 

HD 82286 = FF UMa. A double-lined binary with 
a short 3.275-d orbit and a very precise mass ratio 
of 0.4632±0.0002 is given. Both component's spectral 
lines are rotationally broadened to nearly the same value 
(38 kms^if] and appear variably asymmetric, most likely 
due to cool surface spots. Our photometry shows an am- 
plitude of 0.16 mag with a clear and precise period of 
3.27629±0.00003 d. We attribute it to the brighter primary 
component (Une ratio of 2.9). The spots are the largest con- 
tributor to the radial velocity error budget in our orbital so- 
lution, in particular for the secondary component. Residual 
periods of 1 .945 d and 8.714 d are evident from its radial ve- 
locity. Given the 9.6-mas parallax and the 7.8-mag system 
brightness, we expect two subgiant stars, already noted by 



Henry et al. ( 1995 1 and in our paper I, while Galvez et al. 
( |2007| l favored the odd KOVh-KIIV classification with the 



' In paper I we measured vsini for th e two components of 35 and 
29 kms^^, respectively, but Galvez et al. j2007} quoted us with values 
of 17 and 16 kins~^ which we do not know where they came from. 
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Fig. 9 Seasonal orbital periods for HD 82286 (FF UMa). 
Shown is the period difference in days versus truncated Ju- 
lian date for four independent orbital solutions with respect 
to the period that is obtained from the full data set (i.e. 
3.2751924 d; see Table |5]l. 



less-massive component to be the evolved star. Two sub- 
giants are more likely and would favor the 1.945-d period 
to be the true rotation period of the secondary (suggesting a 
minimum radius of 1.45 R©). The Ha line appears in emis- 
sion and with a strongly variable profile (Fig.[8]l. 

Galvez et al. ( |2007[ ) claimed to have detected an orbital- 
period variation on the order of dP/ P —5 10^*, but Griffin 



(20121 could not confirm this variation from a reanalysis 
of the same data. Due to this situation our STELLA data 
received a special treatment of the systematic errors. We 
followed the procedure outlined by Weber & Strassmeier 



(20111 for Capella. Briefly, we compute synthetic spectra 
for both components and combine them according to their 
real brightness ratio and then extract a radial velocity with 
the same 2d-cross correlation as for the real data. This en- 
ables the removal of the systematic blending effects and 
converges on rms residuals smaller than the precision of 
a single data point. With this, our four consecutive years 
of observations do vaguely hint towards such an orbital- 
period variation but unfortunately remain still inconclusive 
(Fig. [9]). We computed orbits for all four observing sea- 
sons separately as well as for the entire data set and de- 
rived upper and lower limits for dP/dt between 2.76 10^^ 
and 0.2 10^^. In case we take an equal-weight approach for 
the four solutions the most likely value for dP/ P would be 
5.5 10^^, with an upper limit of L2 10^^ and a lower limit 
of 8.5 10^^, i.e. close to zero; altogether not quite consistent 
with above claim by Galvez et al. (2007). Fig. |9] shows the 
orbital periods of the four STELLA observing seasons over 
time. 

HD 82841 = OS Hya. Our photometry shows a period of 
5L54±0.02 d while the SBl orbital solution gives a period 
of 56.7316 d with a slightly eccentric orbit with e=0.087. 
The pseudo-synchronous rotation period would be expected 
to be near 53.9 d. As the photometric period is even shorter 
than this, the star is likely an asynchronous rotator. Our 



V sin i measure of 8 km s^^ converts to a minimum radius of 
8.1 R0 and, together with the temperature and gravity from 
our spectrum synthesis, suggests a K2 giant. We note that 
the metallicity of -0.30±0.03 is significantly subsolar. 

HIP 50072 = EQ Leo. We give an SBl orbit with 
a formally but marginally significant eccentricity of 
0.014±0.0043 and a period of 34.3157 d. An e=0 solution 
has an rms of the residuals worse by salO % . The Ha-line 
profile appears as an asymmetric and variable absorption 
line, indicating chromospheric activity. The Ha-core flux 
shows a period of 34.8 d, very close to the orbital period 
while the photometry gives a 33.36 d as the best-fit period. 
The latter is better defined and we interpret it to be closest 
to the true rotation period of the star. The system is then in 
or close to dynamical equilibrium. 

HD 93915. Our SB2 orbit from STELLA data is in good 
agreement with the Cambridge orbit from Griffin ( |2009| l. 
As our baseline is more than five times longer and better 
sampled, our orbit clearly supersedes the Cambridge orbit. 
The orbital period is near 223 d and e=0.3795±0.0009. The 
expected pseudo-synchronization period is will d but our 
photometry does not show any type of periodic behavior de- 
spite of having a dispersion of up to 0™2 mag in V over 138 
days. 

HD 95188 = XZ LMi. We judge this star to be single. 
The radial velocities possibly show small-amplitude varia- 
tions of ±0.1 km s^^ but no clear period while the Ha-core 
flux indicates a period of 190 d, most-likely spurious. Our 
photometry revealed an uncertain period of 7.00±0.06 d but 
with a small amplitude of just 0.025 mag in V , that we nev- 
ertheless interpret to be the rotation period. 

HD 95559 = GZLeo. The star was found to be a double- 
lined spectroscopic binary by Jeffries et al. ( jl994j who also 
predicted its photometric variability. The orbit was later re- 
fined and the photometric variability discovered by Fekel 
& Henry ( |2000j . Our new circular orbit is very well de- 
termined despite the rotationally broadened line profiles of 
vsmi of 33 kms~^ for both components. In our original 
survey paper (paper I), we noted that the star may be triple- 
lined based on a single coude spectrum. We find no evidence 
for a third component in our new echelle spectra. However, 
the secondary's residual radial velocities of 0.165 kms~^ 
indicate a period of 0.592 d that appears to be an alias of 
the 1.517-d period obtained from our photometry. The true 
rotation period is likely 1.517 d as it also agrees with the 
orbital period of 1 .526 d. Obviously, the system is well syn- 
chronized. 

HD 95724 = YY LMi. Its light variability was discov- 
ered in our paper I and accordingly named YY LMi by 
Kazarovets et al. (2006). The star has a constant radial ve- 
locity from 76 STELLA measurements and is likely single. 
Our APT photometry shows a period of 1 1 .49 d and con- 
firms that YY LMi is a spotted and active star. 

HD 104067. Single star with no obvious periodicity 
down to the detection limit of «100 ms^^ for the given 
vs'mi of 8 kms^^. However, there are two equally signifi- 
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cant peaks in the periodogram from the Ha flux of 34 d and 
102 d, one the alias of the other. We had no photometry for 
this target. 



HD 105575 = QYHya. Szalai et al. ( 2007 1 presented a 



photometric and spectroscopic study of this 0.29-d eclipsing 
binary and derived a complete system's solution from BV 
light curves and radial velocities. These authors also an- 
nounced the discovery of a tertiary component in the cross 
correlation function and listed absolute parameters appro- 
priate for a (K5/M1)/G4 system. We dare to note that they 
(and the referee) had missed our paper I in A&A from 2000 
where we gave radial velocities of the binary plus the ter- 
tiary and determined Ha emission-line fluxes for all three 
components. We then noted that the primary had moderately 
strong and sharp Ca II H& K emission and showed that the 
resonance lines appeared to be asymmetric due to the reso- 
nance lines of a third component. The Ca II emission lines 
even appeared with sharp and red shifted interstellar absorp- 
tion. Because our new STELLA data sample only one fourth 
of the tertiary orbit, we added our earlier paper-1 KPNO 
coude-feed spectra, two HARPS and one FEROS velocity 
from the ESO archive. With these, we extent the time cov- 
erage to 4489 d and obtain an orbital period for the tertiary 
around the eclipsing pair of 4550±52 d and give a prelim- 
inary orbit with e=0.391±0.006 for it. Our close-pair orbit 
is SB2 and circular with a period of 0.2923378 d. Its mass 
ratio is 0.457. The secondary's K amplitude is 158 kms^^ 
and both stars have a wsini of around 40 kms^^; and ac- 
cordingly large are the orbital residuals. The disagreement 
with the elements given by Szalai et al. ( |2007 [) possibly 
comes from the different ways the cross correlation is de- 
termined. Our photometric orbital period is 0.292386 d and 
agrees with the orbital period within its uncertainties. The 
eclipse depth is 0.36 mag in V. No rotational modulation is 
detected from the radial velocities. 

HD 106855 = UV Crv. The catalogue of compos- 
ites of double and multiple stars (Dommanget & Nys 
[199T) identifies HD 106855 as the double star Don 520 
= CCDMJ12174-2104AB. The B component is just 1.3" 
away from the A component at a position angle of 38° 
(both values measured in 1959 by Donner) and appears ap- 
proximately 1.7 mag fainter in V than the A component. 



Worley (1972l gave a more recent measurement of 1.64" 
and 31.6° in 1972.1. The Hipparcos/Tycho photometry lists 
Vt=9.76±0.03 and {B - 1/)t=0.95 for the A component 
(Fabricius & Makarov |2000| l. The fiber-core diameter of our 
STELLA/SES projected on the sky is 3.0" centered on the 
brighter A component. Therefore, the B-component is just 
at the rim of the field of view and could have been within the 
diaphragm from time to time. Our cross-correlation function 
shows a broad, sometimes doubled, main peak and a third 
narrow peak which is sometimes stronger than the main 
peak. 

The orbital period from the main peak, i.e. presumably 
from the Aa component, is 0.67 d with strong aliasing at 
«2 d. There are 1 1 phases where the main peak is split into 



two peaks and possibly represent detections of the Ab com- 
ponent (plotted in Fig. [3]l. However, these velocities have 
unreasonable high scatter (9.9 kms^^) and a comparably 
small X-amplitude (37.7 km s^^). Possibly, the Ab compo- 
nent is part of another spectroscopic binary. Its velocities 
alone would suggest a 0.2-d period. A tentative 0.85-d pe- 
riod is found from the third, narrow peak of the cross corre- 
lation function but we can not claim it conclusively. How- 
ever, it clearly appears sometimes on the one or the other 
side of the main peak, but not always. One possibility is 
that the visual B component is only sometimes within the 
STELLA fiber aperture and is itself a spectroscopic binary. 
It would make HD 106855 a quadruple system. 

HD 108564. Single star with a high average radial ve- 
locity of H-1 1 1.363±0.009 kms^^. A periodogram of the 
radial velocities hint at a period of 9.58 d with a small 
amplitude of 0.1 kms^^ while the Ha flux favors a (for- 
mally even more uncertain) period of 3.63 d. The vsmi of 
14.8 km s^^ would convert to a minimum radius of wl 
if 3.63 d is interpreted as the rotation period. The 9.58-d 
period would result in a minimum radius in disagreement 
with the Hipparcos distance of 28 pc. Our spectrum synthe- 
sis clearly indicates a late K dwarf with rcff=4560±24 K, 
log (7=4.40±0.06 and a significant sub-solar metallicity of 
-0.90±0.03 which implies a radius of «0.7 R©. 

HD 109011 = NO UMa. Griffin (20TF) noted for 
this star "Slow changes in velocities, SB2". We verify 
the system to be a SB2 and find an orbital period of 
1274.7±0.9 d. Our photometry indicates variations with a 
period of 8.4±0.2 d and an amplitude of 0.02 mag in V, 
which we interpret to be due to rotational modulation of 
the primary component. This period is similar but not com- 
pletely in agreement with the values of 8.56 and 9.19 d 



found by Gaidos et al. (20001. But in any case, the system 
contains at least one very asynchronously rotating compo- 
nent because the expected pseudo-synchronization period 
would be «440 d (e=0.507). The orbital solution shows 
radial-velocity residuals of 0.17 and 0.46 kms~^, respec- 
tively, which are larger than expected for the v sin i given. 
A Lomb-Scargle periodogram of these residual velocities 
indicate periods of 639 d for the primary and 246 d for the 
secondary, the former roughly half the orbital period. The 
large Hipparcos parallax and a 1^=8.1 mag system bright- 
ness suggests two dwarfs rather than giants. If interpreted 
due to rotation, the periods of w639 and 246 d would in- 
dicate minimal radii appropriate for supergiants, totally im- 
possible for a distance of 23.7 pc. If these periods are real 
at all, they are possibly connected with surface activity, e.g. 
due to a spot cycle or active longitudes. The adopted ro- 
tation period of 8.4 d and the wsini of 5 kms^^ for the 
primary suggests a minimum radius of 0.83 R©. 

HD 111487 = IM Vir. An SB2 orbit with a pe- 
riod of 1.3086147 d is given. Our photometry confirms 
the eclipsing nature of the system with a period of 
1.3085900±0.000007 d and an eclipse depth of 0.76 mag 
in V. A rms of the residual velocities of the primary 
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Fig. 10 Ha line profile of HD 11 1487 at two orbital 
phases 180° apart. Its secondary component is detectable 
only in Ha where it is seen as an emission line (appearing 
bluewards of the primary's absorption line in the top spec- 
trum, and redwards in the lower spectrum). 

of 0.268 kms~^ is appropriate for rotationally broad- 
ened spectra of t;sini=42 kms^^ but the residuals for the 
secondary-star velocities are anomalously high (1.8 km s^^; 
vsinz=28 kms~^), indicating an asymmetrically spotted 
star However, residual periods of 3.45 d or 1.41 d (one the 
alias of the other) are seen for the primary but none for the 
secondary. Their amplitude is wO.6 kms^^. The secondary 
shows strong Ha emission while the primary appears with 



Ha in absorption (see Fig. 10 1. The photometry reveals a 
period of 1.309 d from data outside of eclipse and have a 
full amplitude of 0. 15 mag in V that we interpret to be from 
rotational modulation synchronized to the orbital motion. 

HD 1 12099. A single-lined eccentric orbit with a pe- 
riod of 23.50447±0.00013 d and an rms of 67 ms^^ is 
obtained. The expected pseudo-synchronization period is 
«9.6 d (e=0.447). A residual period of 4.6 d is seen in a 
Lomb-Scargle periodogram but is not conclusive. Together 
with our Tjsini measure of 4.2 kms^^, it would indicate 
a minimum radius of 0.38 R0, at least in principle agree- 
ment with the Hipparcos parallax. Our spectrum synthesis 
confines this to a KIV classification, in agreement with its 
B~V color. 

HIP 63322 = BD+392587B. B component of a visual 
pair with HD 1 12733 (HIP63317) as the primary. The spec- 
trum appears at almost all wavelengths single lined but 
shows doubled Call H&K lines (Strassmeier et al. |2000| ). 
Galvez et al. (12006)1 presented an SB2 orbit from data taken 
during nine nights in March/ April 2004. They found an ec- 
centricity of 0.3 and noted that the spectrum is best fit with 
a G8V primary and an MOV secondary star An SBl orbit 
with zero eccentricity and P=4. 593843 d is derived from 



STELLA data. Our spectra show the secondary with Ha in 
emission and the primary with Ha in absorption, similar to 
HD 1 1 1487. The Ha-line flux of the primary shows a clear 
periodicity of 2.30 d, exactly half the orbital period within 
the uncertainties, which is simply due to the fact that the 
secondary-star emission is sweeping through the primary- 
star line. The photometry reveals a period of 2.227±0.003 d, 
assumed to be from a double-humped light curve, and we 
adopt twice its value to be the true rotation period of the 
star. 

HD 112859 = BQ CVn. Just recently. Griffin p009] l ob- 
tained an orbit very similar to ours and we refer the reader 
to his discussion on the various stellar properties. Our orbit 
was determined from 93 velocities for the primary and 87 
for the secondary over a time range of 1249 d and achieved 
rms residuals of 0.30 kms^^ and 0.26 kms^^ for the pri- 
mary and secondary, respectively. Our photometry indicates 
a rotational period that is synchronized to the orbital revo- 
lution. 

HIP 63442 = CD CVn. A single-lined orbit with 
P=38.7156 d and e=0 is given. We see many detections 
of a second line system in the cross-correlation functions 
but with such a large scatter (9.9 kms~^) and small ampli- 
tude (0.5 kms^^) that we can not be sure of its reality. Its 
marginally significant period is in agreement with the or- 
bital period of the primary but in phase with all its measure- 
ments. It resembles HD 16884, possibly being a pair of SBl 
binaries in a bound orbit. The residuals from the primary 
orbit have an rms of 0.27 km s~^ and show a periodicity of 
19.14 d, nearly exactly half the orbital period. The Ha-core 
flux varies with a period of 39.2 d and our photometry gives 
a period of 40.3 d, which we interpret to be the true rota- 
tional period of the star. Thus, CD CVn is rotating not quite 
but close to synchronous. 

HD 127068 = HKBoo. We present a new SB2 orbit with 
a period of 15.15011±0.00011 from 98 velocities for the 
primary and 99 for the secondary spanning a time range of 
1257 days. The residual rms is 0.19 kms^^ for the primary 
and 0.36 km s^^ for the secondary. There are periods in the 
residuals of both stars of 18.6 d and 24.3 d for the primary 
and secondary, respectively. The primary's value is in agree- 
ment with our photometric period of 18.187 d in this paper 
and suggests that it is likely the rotational period of the star. 
The secondary's period is maybe spurious as its amplitude 
is also very low, 0.3 km s^^. From photometry in paper I we 
found 17.62 d. In any case, this would make the HK Boo pri- 
mary an asynchronously rotating star Griffin (2009 ) gives a 
discussion of the star's properties and cited a photometric 
period from Hipparcos data of 6.953 d that we can not ver- 
ify. 

HD 136655. In paper I, we had given two velocities of 
-29.9 ± 4.1 and -1-17.9 ± 1.7 kms^^ from a single blue- 
wavelength spectrum for two components. The latter must 
have been an erroneous measurement because it is not ver- 
ified in the present paper Table [3]lists an average STELLA 
velocity of —31.921 ± 0.035 kms^^ from 122 spectra. Its 
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Ca II H& K emission appears very weak, if detectable at all. 
Our photometry shows the star to be constant. 

HD 138157 = OX Ser. Our SBl orbit suggests a very 
small eccentricity of 0.005 1±0.0009 from a total of 148 ve- 
locities. The rms of the residuals is 0.22 kms~^, an excel- 
lent value for the large line broadening of v sin i=4Q km s^^, 
and somewhat better than if a circular orbit is enforced (rms 
of 0.25 kms^^). A peiiodogram from the Ha-line fluxes 
gives the strongest peak at a period of 14.38 d, very close 
to the orbital period. Our photometry shows clear variations 
but suggests a period of 7.18 d, pretty much half the period 
of the orbital revolution. Its phase stability indicates that the 
amplitude is dominated by the ellipticity effect but annual 
variations indicate that part of it is due to spots. We con- 
clude that the rotation period of the star is very near the or- 
bital period and thus synchronized with the orbital motion. 

HIP 77210 = V381 Ser. The star's alternative name is 
G16-9. An initial SB2 orbit with a period of 9.94319 d 
was determined by Goldberg et al. ( 2002| l. Our SB2 orbit 
is among our most precise with an rms of the residuals of 
the primary of 52 ms^^ from 96 velocities (but four times 
larger rms for the secondary). Our photometry shows vari- 
ations with a period of 13.7±0. 1 d that we interpret to be 
the rotation period of the primary, thus, significantly asyn- 
chronous with respect to the orbital period. The secondary's 
velocity residuals show a systematic rms with a period of 
4.79 d, a little less than half the orbital period. The orbital 
eccentricity of 0.06032±0.00024 suggests a formal pseudo- 
synchronization period of sa9.7 d, which is close to twice 
the secondary's residual period. We conclude that 9.58 d is 
likely the true rotation period of the secondary and thus syn- 
chronized to the orbital revolution. 

HD 142680 = V383 Ser. The photometric period from 
data in 1999 in paper I was 33.52 d. Our new photometry 
principally verifies this period and gives 33.4±1.5 d. The 
orbital eccentricity is 0.3158 from our new STELLA SB2 
orbit and is in agreement with the eccentricity of 0.3139 
from Griffin's (2009) SBl orbit. The pseudo-synchronous 
rotation period would be around 15.3 d, which is suspi- 
ciously close to 1/2 of the photometric period. However, the 
photometric period can not be halved if it is due to spots 
rotating in and out of view and we conclude that the sys- 
tem is rotating super synchronously. The systemic velocity 
is relatively high with —83.83 kms~^, making this system 
a high- velocity binary. 

HD 143937 = V1044 Sco. This is the A component of 
the visual binary GJ 9539. The B component is 10.3" away 
and approximately 2.7 mag fainter than component A in 
the Hipparcos bandpass. It is therefore included in our pho- 
tometry with a 30" diaphragm but well separated for the 
STELLA spectroscopy. Cutispoto et al. ( 1999) found the A 
component to be an SB2 binary (it is listed as an eclipsing 
binary in MaUcov et al. |2006| based on a note in Cutispoto et 
al. |1999| . We verify its eclipsing nature and find an orbital 
period from photometry of 0.9 1479±0. 00003 d, in excellent 
agreement with the radial-velocity orbit of 0.91483076 d. 



The eclipse depth is 0.41 mag in V while the spot wave 
shows an amplitude of 0.10 mag. The rms of our orbital 
residuals for both stars are comparably high, 0.28 kms^^ 
and 1.26 kms~^, but appropriate for the vs'mi of 55 and 
45 kms~^ for primary and secondary, respectively, jsucvh 
large residuals are related to the short orbital period of 0.9 d 
and the corresponding degree of line blending. 

HD 147866 = V894 Her. We present an SBl orbit with 
a period of 114.050 d and an eccentricity of 0.5141. The 
K amplitude is just 6.71 km s~^ and the mass function ac- 
cordingly small. Our photometry reveals a likely rotation 
period of 80.3±2.2 d, while pseudo synchronization would 
be expected at a period of ss40 d. It appears that the system 
had not had enough time yet to obtain orbital and rotational 
equilibrium. With our vsini measure of 5±1 kms^^ the 
minimum radius is 7.93 Rq. The Hipparcos parallax places 
the star at a distance of 193 pc which implies an absolute 
magnitude of H-1.67 mag (for V=8.1 mag). Our spectrum 
synthesis gave Toff =4580 K and log 5=2.7 with a significant 
sub-solar metallicity of -0.24±0.02. The star is most likely 
a KlIII giant. 

HD 150202 = GI Dra. We confirm Griffin's (|2009} 
orbital period with a value of 68.47160±0.00097 d. The 
STELLA orbit's rms of the residuals is just 83 m s^^ from 
179 velocities over 1 1 12 d. Our spectrum synthesis yielded 
Toff=5010±25 K and a log.g=3.06±0.07, suggesting more 
a G8III star than a KOIII as put forward by Griffin ( p009] l. 
The periodogram from our photometry shows the strongest 
peak at a period of 36.4±0.5 d. This is likely just one-half of 
the true rotation period (72.8 d) and is then in rough agree- 
ment with the period of 76.70 d from Hipparcos photometry. 
Thus, the star is an asynchronous rotator with Prot > ^orb- 

HD 153525 = V1089 Her. From 196 STELLA spectra 
we obtain a mean radial velocity of -6.844±0.039 kms~^. 
The star is single but at least moderately active judging 
from its Call H&K emission-line strength (c/o paper I). 
The new photometry shows variations with a period of 
15.4±0.1 d that we tentatively interpret to be the rota- 
tional period. A distance of just 17.5 pc and our values for 
reff=4775±23 K and log5=4.47±0.05 suggest a K dwarf 
with a V sin i <3 km s~^, in agreement with our upper limit 
of 3 kms^^. 

HD 153751 = e UMi. The system is a single lined spec- 
troscopic binary and has an orbit determination dating back 
now 60 yrs (Clime nhaga et al. |1951j ). This is certainly be- 
cause of the high declination of +S2° despite that Climen- 
haga et al. already noted the Ca II H&K emission and thus 
put the system out for further observation. Its eclipsing na- 
ture was discovered by Paul Guthnick in 1946 and later ver- 



ified by Hinderer ( 19571 who concluded on a pair of gGl 
and dA8-dF0 stars for the primary and secondary, respec- 
tively. Today's accepted classification is G5III for the pri- 



mary (e.g. Keenan & York T988JI. Our spectrum synthesis 
gives Teff=5215±47 K, logg=3.21±0.08 and a metalhcity 
of -0.25±0.04, i.e. significantly subsolar but altogether in 
agreement with a 05111 classification. Our new SBl orbit 
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has a period of 39.48042±0.00012 d and zero eccentricity. 
Its rms of 121 ms^^ from 220 velocities spanning 2144 
days is mostly due to the moderately large rotational ve- 
locity of 25.6 kms^^ and only less from the jitter from 
starspots because no convincing rotation period is found 
from the residuals. 

HD 155802. In paper I, we had given veloci- 
ties of -33.5±1.2 and -33.3±0.5 kms^^ from a blue- 
wavelength and a red-wavelength spectrum, respectively, 
and H-24.7±2.0 km s~^ for another component from the one 
blue-wavelength spectrum. The latter must have been an er- 
roneous measurement because it is not seen in the spectra in 
the present paper The average velocity from 1 14 STELLA 
spectra is -32.71 1 ±0.033 kms^^. Its rms is practically the 
precision of the STELLA/SES system and thus the star is 
constant down to our detection limit of «30 ms^^. How- 
ever, its chromospheric activity is verified by the presence 
of weak Ca II H&K emission lines. 

HD 184591. The star is single according to our 155 
STELLA velocities. Its mean is -37.635±0.044 kms^^. 
The photometry indicates just a weakly determined period 
of P=48.9±2.2 d with an amplitude of below 0.01 mag 
in y. Our spectrum synthesis resulted in Toff=4800±40 K, 
log5=2.82±0.07 and a metallicity of -0.18±0.03. If the 
48.9-d period is real and the true rotational period of the 
star, the v sini of 4.7±0.8 kms^^ confines the radius to be 
at least 4.5 R©, in principal agreement with above gravity 
and the Hipparcos parallax. 



HD 190642 = V4429 Sgr. Strohmeier (TwT) discov- 
ered this star on his Bamberg plates to be variable with an 
amplitude of 0.25 mag. No other information was given. 
But based on the Hipparcos satellite data, Kazarovets et 
al. (2006 1 assigned a variable star designation. We find 



the star to be an eclipsing system and detect a photo- 
metric period of 24.8339±0.0009 d after excluding the 
eclipse points which, if doubled, is very close to the or- 
bital period from radial velocities (49.5969±0.0012 d, 
e=0.0126±0.0011). The eclipse depth is 0.43 mag. The rms 
(0.445 kms^^) of the velocity residuals is larger than ex- 
pected for V sin «= 19.6 kms^^ and a periodogram analysis 
shows a convincing residual periodicity of 24.72±0.02 d 
with an amplitude of 1.2 kms^^. The Ha-core flux shows 
also a significant period at 24.9 d. All of these periods must 
be doubled to obtain the true rotation period and then sug- 
gest that the primary's rotation is synchronized to the or- 
bital motion, or at least nearly so. Our spectrum synthesis 
gives values of Tcff=4760±63 K and log.g=2.96±0.11 with 
a significant sub-solar metallicity of -0.47±0.05. The Hip- 
parcos parallax suggests an absolute system magnitude of 
H-1.34 mag (V^=8.08 mag, distance of 223 pc, no extinction). 
Altogether these values are compatible with a KlIII-IV gi- 
ant. Knowing i «90° gives directly a radius of 19.6 Rq 
and a rotation period of «50 d, in good agreement with our 
various observations. Assuming a mass of two solar masses 
for the primary, the mass function of 0.330 demands a sec- 



ondary mass of 1.22 M0; it suggests a late F main-sequence 
star, maybe F6-8. 

HD 197913 = OR Del. This is the A component of the 
visual binary ADS 14270. Griffin ( 2005[ l gave an overview 
of its observational history and presented a full SB2 or- 
bit in very good agreement with our present determina- 
tion. Our rms residuals are very low though, 100 ms~^ 
for the primary and 161 ms~^ for the secondary. In pa- 
per I, we presented double-lined Li spectra for the A com- 
ponent and a single-lined spectrum for the B component 
(erroneously identified with lower-case a and b in that pa- 
per; pace Griffin). The fainter B component appeared with 
somewhat stronger Call H&K emission as the (blended) 
A components. The 6.56-d photometric period that we de- 
rived from our APT time series from 1999 was assigned to 
HD 197913, without an A nor a B specification. This was 
because our APT diaphragm was 30" and the AB separa- 
tion roughly 11". Consequently, it prompted Kazarovetz et 
al. 



p006 



to assign the variable star designation OR Del 
to the AB system but Simbad resolves it for the A com- 
ponent of which, however, it does not give the HD num- 
ber. Variability from the B component can not be excluded 
but it is indeed more likely that the photometric variability 
stems from the 0.8-mag brighter (combined) A components. 
Our current photometry has a best-fit period of 6.56±0.07 d. 
The AB period is suggested by Griffin (2005 ) to be around 
"thousand years". Our new STELLA orbit gives a period of 
1778. 5±1. 3 d for the A component from basically a single 
revolution (time span of data is 1848 d) but still better con- 
fined than the published value of 1790±13 d. 

HD 199967. We find it to be an SB2 in a 6.100213- 
day orbit with a very low but significant eccentricity of 
0.0024±0.0002. An e=0 solution has an rms of the residuals 
worse by more than 10 % for each component. No photo- 
metric period is detectable but such a short orbital period 
makes it likely a synchronized system. Our v sin i measures 
of 11 and 9 kms^^ for the primary and the secondary, re- 
spectively, then lead to minimum radii of 1.3 and 1.1 R0. 
The effective line ratio is 3.2 to 1. The Hipparcos parallax 
gives a distance of 74 pc and an absolute system magnitude 
of H-3.67 mag, basically excluding two main-sequence stars. 
The system's B — V is 0.58 mag and, with a mass ratio of 
0.776, does suggest two early G-type dwarf-subgiants. 

HD 202109 = C Cyg. The primary is a mild bar- 
ium star with a white-dwarf companion. Griffin & Keenan 



( 1992 1 discussed its spectral classification in detail and pre- 
sented an orbit with a period of 6489±31 d, e=0.22±0.03 
and an amplitude of isr=3.31±0.12 kms^^. More re- 
cently, Eaton & Williams ( |2007| hsted a mean velocity of 
19.55±0.25 kms^^ from 74 observations but did not give 
the individual values. Our own radial velocities show a clear 
non-linear trend with a full range of 2.8 kms^^ during the 
3.5 yrs of observation. We present a new orbit from the com- 
bined values from CORAVEL, the Cambridge spectrometer, 
STELLA and the literature cited in Griffin & Keenan ( |1992| l 
in Fig.[2]and list it in Table[4] It covers a time range of 1 14.8 
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years! We shifted the Griffin & Keenan (1992i points by 
+0.50 kms"^ to bring them from the lAU to the STELLA 
system. The STELLA data in the panel are the ones between 
phase 0. 1-0.3 and an average rms of 42 m s~^. The refined 
period is now 6446±14 d and e=0.244±0.016, in fair agree- 
ment with the Griffin & Keenan orbit. 



HD 218739. Visual component to HD 218738 = KZ 
And, itself a cataloged RS CVn binary 15.6" away. 
HD 218739 appears with a constant radial velocity of 
-5.201±0.040kms-i. 



HD 226099. Duflot et al. ( fT995] l listed six velocities be- 
tween +22 and +4 kms~^ which were the first measure- 
ments that showed the variable velocity. In our paper I, we 
gave velocities for both stellar components from a single 
blue-wavelength spectrum but did not find significant Ca II 
H&K emission lines. In this paper, we present its first or- 
bit (SB2) with a period of 18.78209 d and an eccentricity 
of 0.3096 (Fig. [3]). The two components are of compara- 
ble mass with a ratio of 0.897. The v sinz for the primary 
component is below our resolving limit and estimated to 
be 2 kms^^, but twice as larger for the cool secondary. 
The rms residuals for the two components are 85 m s~^ and 
182 ms~^, possible due to the sharp lines. The orbital ec- 
centricity suggests a formal pseudo-synchronization period 
of w 1 1 .6 d but no we have no photometry for this star. How- 
ever, the radial-velocity residuals of the secondary indicate 
a well-defined period of 3.131±0.008 d with a full ampli- 
tude of up to 1.6 kms"^, and no power for periods longer 
than ss4.5 d. We conclude that the 3.1 days is the rotation 
period of the secondary and that it rotates significantly asyn- 
chronous. 



HD 237944. The system has a triple-lined spectrum due 
to Aab and the visual B component and, according to Otero 
& Dubovsky ( 2004| l, is also an eclipsing system (ab compo- 
nents). The two stronger lines from component A make up 
the ab SB2 in our Table |5] The third, visual B or spectral c 
component's average radial velocity is -12.25±1.3 km s^^, 
different by 3.0 km s^^ from the systemic velocity of ab, in 
accordance with an expected orbital period of «2500 years 
noted by Griffin ( |2009) l. Our rotational velocities are from a 
single (well-exposed) spectrum and gave v sin i of Aa, Ab, 
and B of 1 1 km s~^, 9 kms^^, and 11 kms^^, respectively, 
with errors of at least ±2 kms~^. Note that the SB2 orbital 
period is just 5.5 d with a very small but significant eccen- 
tricity. Our photometry confirms the eclipses of the ab pair 
(Fig.|4]i and also shows a sinusoidal light variation with the 
orbital period indicating that a spotted star is rotating syn- 
chronously with the orbit. Therefore, and assuming i = 90°, 
radii are 1.19 Rq and 0.98 R© and masses are 1 .042 M© and 
1.027 Mq, respectively. For our Li abundances in TablejT] 
we assume effective temperatures based on Gray's [2005] ta- 
ble B.l and Griffin's (2009 1 individual absolute magnitudes. 



6 Analysis 

6.1 Rotational synchronization 

Before we try to derive a rotation-activity relation, we 
inspect our binary sample in the rotational-period versus 
orbital-period frame and identify its distribution of syn- 
chronicity, or the deviation thereof. Asynchronously rotat- 
ing binaries, or pseudo-synchronously rotating binaries in 
eccentric orbits, could experience enhanced activity due 
to stronger shearing forces in their convective envelopes 
(Schrijver & Zwaan 1991 1. This may possibly be coupled 



with enhanced deep convective mixing as in Hertzsprung- 
gap giants (e.g. Bohm-Vitense |20C)4| . It may thus be reward- 
ing to visualize the Prot-^oib plane. 

Fig. [TT| shows the distribution of synchronism in our ac- 
tive binary sample. The full sample comprises of the entries 
in the "Chromospherically Active Binary Star (CABS)" cat- 



alog (Strassmeier et al. [19881 [T993| Eker et al. [2008) 1, ex- 
panded by the stars from the present paper. Of the approx- 
imately 410 CABS binaries only 230 are known with both 
a rotational period and an orbital period (150 have no rota- 
tional, 30 have no orbital period). We note that care must 
be taken when using CABS-III because, e.g., the photo- 
metric periods listed are mixed values of rotational periods 
and Porb/2-periods due to the ellipticity effect or are sim- 
ply orbital periods in case the system is eclipsing. Besides, 
some of the photometric periods are derived from v sin i 
measurements and radius estimates (e.g. for HD 37171 of 
^rot=300 d) and are accordingly uncertain. From Fig. 11 
we see that the bulk of active binaries is synchronized and 
in circular orbits but from the 47 systems that have ro- 
tational periods considerably shorter than the orbital pe- 
riod 41 are with e > and six with e=0 (HD7205, 
HD16884ab, see Figs. [2[ and [4| HD 72688 HD 90385, 
HD 108102, and 93 Leo). These six are definitely sub- 
synchronous rotators and all have orbital periods less than 
approximately 100 d. Of the 41 systems with e > 0, only 
seven (HD 54563, HD 118981, and HD 133822=HSLup if 
its photometric period is exactly doubled, plus the four sys- 
tems already listed in Hall 1986) are rotating at or near the 
expected pseudo-synchronized value according to the crite- 
ria of Hut (1981), leaving 34 more asynchronous rotators 
in the sample. Among these, the two most extreme cases 
are HD 109011 (F,.ot=8.34 d, Forb=1275 d, e=0.508) and 
HD 197913 (Prot=6.62 d, Porb=1790 d, e=0.265), rotating 
54 times and 190 times faster than pseudo synchronization, 
respectively. If real - see the individual note for HD 197913 
- these systems either had significant mass transfer in the 
past or are so young that magnetic braking was not effective 
yet or they experienced a yet unconsidered spin-up mecha- 
nism in the past, e.g. the engulfing of Jupiter-sized planets 



or brown dwarfs (see e.g. Livio & Soker 2002 Carney et al. 
[2003 [ l. We will come back to this thought later in the section 
on lithium in active stars (Sect. [63] i. 

Stars where the rotational period is longer than the or- 
bital period are in any case asynchronous rotators, indepen- 
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Fig. 11 Rotational period versus orbital period for the bi- 
nary stars from Table 1 (big dots) and other RS-CVn and 
BY-Dra binaries from the CABS catalog (small dots). Open 
symbols indicate circular orbits, a filled symbol indicates an 
eccentricity larger than zero. The lines depict the (pseudo)- 
synchronization periods for various values of the orbital ec- 
centricity, e, according to Hut (1981 ). Our sample mainly 
consists of synchronized binaries with circular orbits but 
26% are asynchronous rotators. See text. 



dent of their orbital eccentricity. The early discoveries were 
summarized in the paper by Hall ( 1986"), e.g. A And and 
39 Cet=AY Cet. Three more systems stick out suspiciously 
above the 45 -degree line in our Fig. 



11 though. HD 31738 



is a double-lined short-period binary with a recently revised 
orbital period of just 0.45 d (and e=0; Griftin |2009) . Its pho- 
tometric period of 4.55 d was derived from an earlier APT 
campaign (Strassmeier et al. |1989| l which, in the light of 
the new orbital period, was likely an alias of the true or- 
bital period and needs revision. HD 193216 (a quintuple vi- 
sual system) is a new entry in CABS -III and its A compo- 
nent a single-lined binary with Porb=418.77 d and e=0.08 
according to Griffin ( |2002[ ). In CABS-III it is erroneously 
listed with the orbital period of 3.62 d of its visual B com- 
ponent, another SB according to Griffin (2002). The rota- 



tion period given by Wright et al. ( 2004[ l was determined 
for the brighter A-component indirectly from two measures 
of its Ca II H&K 5-index and an empirical rotation-activity 
relation, thus, lacks a real observation. The third system is 
HD 1818 09 (P rot=59.85 d, Porb=13.045 d, e=0.04; Fekel 
& Henry 2005 [ ) and is a well-studied asynchronously rotat- 
ing KO giant/subgiant. In total there are 21 rotators with 

-Prot 



> Porb in Fig. 11 



The more significant ones are 
HD 217188 (KOIII, Porb=47 d, e=0.47), HD 170829 (G8IV, 
26 d, 0.18), HD 33363 (KOIII, 21 d, 0.07), HD 71071 (KlIV, 
16 d, 0.13), and the four stai's HD 82841 (KlIII, 56 d, 0.09), 
HD 127068 (G5-8IV, 15 d, 0), HIP 77210 (K2V, 10 d, 0.06) 
and HD 142680 (K0-2V/K7V, 24 d, 0.31) from this paper 
(see Tables]?] and |5]). Despite that these systems are a mix- 
ture of giants and dwarfs, all of them have orbital periods 
shorter than 60 d with HIP 77210 as short as 9.94 d (and 



10" 



10' - 



10' ^ 



10" 



iQ-' 



O Single, main sequence 

O Single, evolved 

o Binary, main sequence 

O Binary, evolved 




mm ° „ 



3500 4000 4500 5000 5500 6000 6500 7000 

Te„(K) 



Fig. 12 The relation between rotation period and effective 
temperature for active single stars (open symbols) and for 
binaries (filled symbols). Our program stars are now mixed 
with other stars and are not highlighted. The two symbol 
sizes indicate a star's luminosity class. Larger symbols are 
for evolved stars of class III and IV, smaller symbols are for 
main sequence stars of class V. The short dashed line is the 
fit to the I-sequence dwarf stars from Bohm-Vitense ( |2007| l 
while the long dashed line is our fit to the weakly active (sin- 
gle) giant stars. The two full lines indicate the fits to 
the single active stars (upper line) and the evolved binaries 
(lower line). See text. 



Prot=13.52 d). This subclass does not quite follow Mid- 



dlekoop & Zwaan's ( 1981| l argument that giants in binaries 
with periods less than 120-200 days will be synchronized, 
nor Boffin et al.'s ( |1993| ) suggestion that they should be gen- 
erally circularized if the orbital period is shorter than sa50 d. 
Our active binaries with Piot > Porb would need an ex- 
tra rotational acceleration in order to get synchronized and, 
at the same time, need to redistribute orbital momentum in 
order to get circularized (but see Verbunt & Phinney |1995| 
for the interaction of tidal circularization and orbital eccen- 
tricity). Note that the very old (metal poor) binaries with 
dwarf components in the sample of Latham et al. ( 2002[ l 
appeared to reach circularization at a transition period near 
20 d. As rotational synchronization is predicted to occur be- 
fore orbital circularization (Hut 1981), all of above systems 
should have e > 0, which is actually the case, except for 
HD 127068 (this paper) and 39 Cet. To summarize, our full 
sample includes w26% asynchronous rotators (61 stars) of 
which about half are rotating slower than in orbital equilib- 
rium and can not be explained by being too young to be syn- 
chronized yet or by uncertainties in the measured quantities. 
We conclude that active binaries must have gone through an 
extra spin-down besides tidal dissipation and suggest this to 
be most likely due to a magnetically channeled wind with its 
subsequent braking torque on the star as described already, 
e.g., by Mestel ( |1968| l. 
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6.2 The dependency of rotation on effective 
temperature 

Most stars in this paper are within the temperature range 
4000-6000 K, i.e. roughly MO to GO. Such stars undergo 
the more magnetic braking the stronger their initial mag- 



netic field (c/o Gray 2005 and the abundant literature cited 
therein). To first order, this should be proportional to the 
effective temperature. Our first step is thus to extract the re- 
lation between rotation period and effective temperature, in 
particular for the active binary sample where we could ex- 
pect a qualitatively different behavior due to the spin-orbit 
coupling. A similar attempt for single main-sequence G and 
K stars led Bohm-Vitense p007| l to find a strong and sys- 
tematic increase of Prot with decreasing Tcff for "active 
stars" but at most a very weak increase for the "inactive 
stars" of form T^'^^. Its cause is not so clear but ascribed 
to an angular momentum transfer from the surface to the 
interior due to some deep mixing. 

Fig. 12 shows the distribution of Piot versus Tcff in our 



full active-star sample. Also shown are the single dwarf 
stars used by Bohm-Vitense ('2007'), drawn from the Saar 
& Brandenburg ( 1999) sample, and the young solar analogs 
by Gaidos et al. (2000). We also add the single stars with 
H&K emission from our paper I where we had a photo- 
metric period determined as well as the single giants from 



Strassmeier et al. ( |1994 1. From the latter sample of 12 single 
giants, only two had a photometric (or S-index) period that 
we interpreted to be the rotation period (S CrB, HD 160952), 
the others were based on v sin i measures and an estimate of 
the radius from the revised Hipparcos parallax. Note that we 
dropped HD 181943 from the sample (single KIV star seen 
almost pole on) because the photometry likely shows a spot- 
variability time scale rather than the rotational period (e.g. 
Fekel & Henry [T995] ). 

In Fig. 12 the four subgroups - single stars on the main 



sequence, evolved single stars, binary components on the 
main sequence, and evolved binary components - appear 
well separated from each other, as expected due to the dif- 
ferent radii and the tendency towards synchronization in bi- 
naries. The longest rotation periods are found for giants, 
the shortest periods for dwarfs, while the binary compo- 
nents are generally more rapidly rotating than single stars. 
For stars with T^tf <4000 K, one has only dwarf stars to 
work with because M giants exhibit mostly pulsation pe- 
riods rather than rotational periods. However, the strong in- 
crease of Prot with decreasing Tcff is obvious in our sample. 
Two simple functional fits to the single stars and the bina- 
ries, respectively, indicate that the steep increase is not only 
seen in single stars but also in binary components. Their 
slopes are similar or possibly even identical, as expected if 
there is a common mechanism. We find 



Prot = (27.9 ±3.2) T^;7.2±0.9 



■-cff 



Pr, 



(27.1 ±3.5) 



-7.0±0.9 



eff 



for singles, 
for binaries. 
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This is to be compared with the active (A) sequence stars in 
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Fig. 13 Average Ha-core flux as a function of effective 
temperature. Our binary sample is shown as large grey dots 
while the single stars are shown as small grey dots. The sin- 
gle stars from our paper I are shown as open circles. The 



single-star sample from Cincunegui et al. (2007i is shown 
as crosses while the stars in the Lopez-Santiago et al. (2010 1 
catalog are shown as pluses. The line is the minimum-flux 



envelope derived by Cincunegui et al. ( 2007 
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Fig. 14 Ha flux measurements of our sample with two 
different central bandwidths (1 A and 1.5 A). The short line 
is a linear regression fit, the other line is just a 45-deg line 
to guide the eye. Both slopes are identical but the offset is 
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A separate fit to the inactive, or only very weakly active, 
single giants and subgiants with Prot >100 d yields a slope 
of T^'^^, very similar to the inactive (I) sequence fit for 
dwarf stars in the Bohm-Vitense (| 2007| l sample, just offset 
towards longer periods. Note that we plot her fit for inactive 
dwarfs as a dashed line in Fig.[T2]for comparison. Our sam- 
ple also shows a void of single stars for a period range of 
W15-30 d for Toff >5250 K, although not overly striking. 
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6.3 The Ha-line core as an activity indicator 

The Ha-line profile is influenced by a large range of at- 
mospheric optical thicknesses, Doppler widths, and photo- 
spheric line strengths, thereby stretching its formation from 
the photosphere and the chromosphere to the transition re- 
gion with the corona. In normal FGK stars. Ha is domi- 
nated by photoionization. Collisions with free electrons are 
a significant source for the cooler stars and the line becomes 
particularly sensitive for mass motions and all sorts of (non- 
thermal) velocity fields (e.g.. Cram & MuUan 1985) 1. More- 
over, circumstellar hydrogen could add an extra component 
to the net Ha absorption profile. For the average quiet Sun 
the core of the line forms at about 1500 km and is thus 
chromospheric in origin. However, extracting just the non- 
thermal (=activity) contribution from the line core is not 
unique and neither is its interpretation. Particular contribu- 
tors are chromospheric inhomogeneities like bright plages, 
chromospheric network variations, localized in- and out- 
flows associated with flares (nano, micro and macro ver- 



sions) and, of course, global winds. Cincunegui et al. (2007 1 
pointed out that despite they had found a strong relationship 
between the mean fluxes in Ha and Ca II H&K the relation 
can break down for individual stars. They suggested that the 
relationship is the product of the dependence on stellar color 
rather than on similar activity phenomena. 

In Fig. 13 we plot the Ha-core flux, obtained as de- 



scribed in Sect. 4.6 versus effective temperature. Our sam- 
ple is made up of the 7 1 components of active binaries and 
18 single stars from Table [T| with temperatures obtained 
from spectrum synthesis. These entries are complemented 
by the stars from paper I, the 109 southern stars in the study 
by Cincunegui et al. ( |2007| l, and the averages for 86 stars 
from the survey of Lopez- Santiago et al. ( 2010| l. Unfortu- 
nately, we have no stars in common with Cincunegui et al. 
and Lopez-Santiago et al. which makes an empirical deter- 
mination of the zero-point shift between the samples merely 
impossible. The same is true for the observations of F8-K5 
dwarf stars by Pasquini & Pallavicini ( |1991[ ), also done in 
the southern hemisphere, while our survey was done in the 
northern hemisphere. Cincunegui et al. (|2007|) were able to 
connect their zero point with that in Pasquini & Pallavicini 
( |1991| ) for the 24 stars in common and found excellent 
agreement, suggesting fluxes good to within 20%. Never- 
theless, we did not add the Pasquini & Pallavicini ( |1991| l 
sample in our Fig. 13 because it adds additional uncertainty 
due to them using a mix ofV — R and B — V colors. A zero- 
point difference between our fluxes and those of Cincunegui 
et al. ( |2007) is obvious for the single stars in the sense that 
our fluxes appear lower by approximately 50%. We note that 



Cincunegui et al. ( |2007 1 and Pasquini & Pallavicini ( 1991| l 
measured Ha-core fluxes from a bandwidth of 1.5 A and 
1.7 A, respectively, while we used a 1.0-A bandwidth. It af- 
fects how much basal flux is being included in the measure- 
ment, which could be significantly different for the cooler 
stars as well as for the subgiants and giants. To quantify this, 
we have remeasured our spectra with a 1.5-A bandwidth and 



compared them to the 1 .0-A fluxes (Fig. 14 1. We found a ba- 
sically constant offset of 2±2 10^ ergcm^s^^A^^ (where 
the error is the rms of the linear fit) and no dependency 
on color nor on gravity. A correction of this amount was 
then applied to the Cincunegui et al. fluxes before plotting 
in Fig. 



13 



Quite different is the case for the stars observed by 
Lopez-Santiago et al. (2010| ). Thek fluxes were obtained af- 
ter the subtraction of bona-fide inactive standard-star spec- 
tra and likely represent only the active part of the chromo- 
sphere. This is markedly different to all other fluxes which 
contain the photospheric flux but also the basal chromo- 
spheric flux. We converted the B — Vsin their tables to an 



effective temperature with the calibration of Flower ( 1996 1 



and plot them as small pluses in Fig. 13 for comparison. The 
Lopez-Santiago et al. fluxes per temperature bin are on av- 
erage an order of magnitude lower than the other measures, 
as expected. 
Fig. 



13 also shows that the binary components follow 



basically the same trend as single stars but exhibit twice as 
much scatter as the single stars. This is expected because 
the temperature errors for the binary components are typi- 
cally also twice as large, in particular for the SB2s. The line 



in Fig. 13 is the minimum Ha flux for single stars found 



by Cincunegui et al. (2007i from a quadratic polynomial. 
It essentially represents the photospheric contribution. Our 
measurements have a number of stars with fluxes that fall 
below this line. These are almost exclusively the weak-lined 
secondary components in SB2s with comparably uncertain 
flux measurements. 

6.4 An activity-rotation relation for magnetically 
active binaries 

Would a binary component with an equal rotational period 
of a single star be more active just because it is a compo- 
nent in a tidally interacting binary? Or is rapid rotation in 
binaries just an unique evolutionary situation because they 
are speed up to stellar rotations that otherwise would never 
be reached during regular single-star evolution? Tidal cou- 
pling, expressed as the ratio of stellar to Roche radii, was 
ruled out as a direct factor for increased activity in binaries 



by Basri ( 1987 1 despite that the evolved binary components 



in his sample were always more active for a given period 
than a single star (Basri et al. |1985[ ). Yet other authors had 
claimed such a dependency (e.g. Glebocki & Stawikowski 
[1988) 1. More recently, DaU et al. ( [2005| l tried to find a rela- 
tion between binarity, magnetic activity, and chemical sur- 
face abundances by comparing observations of two selected 
cool stars in great detail. What they found was a correlation 
between the bisector inverse slope and the activity index 
log Rhk, which both vary in phase with tiny (few m s~^) ra- 



dial velocity variations. In that paper, Dall et al. ( 2005 1 spec- 



ulated towards the existence of an unseen very-low mass 
star or even planetary companion as the cause of the low- 
amplitude velocity variations, but in their follow-up paper 
of the target in question, Dall et al. ( |2010) ) unambiguously 
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Fig. 15 a Average Ha-core flux versus rotation period, b The same flux but normalized to the bolometric luminosity 
(TT|fjj. The grey scale dots are the singles and binary components from this paper. The other symbols are single stars from 
sources as indicated. The lines are power-law fits as described in the text. 



verified these to be due to stellar oscillations with an w50- 
min period. We conclude that a straightforward comparison 
of two well-selected stars, even with (nearly) equal rotation 
periods, may not be conclusive. 

From statistical attempts, weak activity correlations 
with rotation for binaries were reported in the literature. 
Simon & Fekel (T987') found a marginally decreasing UV 
flux with increasing period for evolved binaries. Strassmeier 
et al. ( |1990| l claimed a slight trend for Ca II H&K, while 
Dempsey et al. ('1993') saw a similar trend for the Ca II in- 
frared triplet. Strassmeier et al. (1994 ) determined a relation 
for Ca II H&K flux in giants versus rotational velocity and 
temperature and emphasized the large range of fluxes for 
a given rotation rate. Also, the overall consensus was that 
the flux for evolved binaries with periods longer than 10 d 
is fairly constant and only increases significantly for stars 
with periods below wlO d. 



In Fig. 15 1, we employ the Ha-core flux, J^Ha. as an ac- 



tivity measure and plot our program stars in the log-normal 
activity-rotation plane. The sample is expanded by those 
stars from paper I that had a photometrically determined ro- 
tation period. The figure shows that the Ha-core flux also 
levels out with longer periods; for single stars near a period 
of 10 d, for binaries more near periods of 20 d, compara- 
ble to what has been seen from other activity indicators. It 
also shows that binary components appear generally more 
active than single stars for a given period bin. Notice that in 



Fig. 15 only those SB2s are plotted where the photometric 
rotation period can be assigned unambiguously to the pri- 
mary star. The secondary stars were not included in the fit 
for the spectroscopic binaries either because their rotation 
periods are assumed from v sin i and a radius instead of ob- 
served, as for the primaries. The lines in Fig. [15^ are simple 
power-law fits of form F = a . 



Fig. 15 3 shows the same flux as in panel a but normal- 
ized to the bolometric luminosity crT^f^, 



cff 



(8) 



Effective temperatures for the majority of stars were de- 
termined spectroscopically from our SES spectra (see Ta- 
ble |7]) or, for the single stars from paper I, from the color- 
temperature calibration from Flower ( | 1996) 1. We also in- 
clude in the single-star fit the targets from Cincunegui et 



al. ( 2007 1 whenever a rotation period from photometry was 



known and an effective temperature given. These data are 
plotted as crosses. Note that the CABS stars have no Ha- 
flux determination as an ensemble but just individual mea- 
surements collected from the literature with a variety of cal- 
ibration issues and are not included here. Again, the lines 
are power-law fits of the form R ~ a P^. The results from 
Fig. 15 1 and Fig.[T5|^ are 

3-0.30±0.05 



j-H„ =(4.22+[;:«io+«)p, 

^Ha = (2.99;°;^3 10+6) p^. 



+0.67 
-0.22 



10- 



i?Ha = (9.72 
i?Ha = (5.92t°;|^ 10 



rot 

0.24±0.04 

■ot 

0.24±0.03 
ot 

0.14±0.03 
ot 



binaries, (9) 
singles, (10) 
binaries, (11) 
singles. (12) 



The rms values for the fits in Eq. |9]-12 are, from top 
to bottom, 0.920eH-6, 0.522eH-6, 1.353e-5, and 1.105e-5, re- 
spectively. We see that the fluxes are only weakly depen- 
dent on the rotation period, in particular when compared to 
the p-i 33 dependency for inactive single stars from Call 
H&K (Bohm-Vitense 2007 ). This is partly because the mag- 
netic flux seen in the Ha core is dispersed over a larger 
range of optical depths at which a variety of (local) veloci- 
ties prevail. However, our fit demonstrates the existence of 
a period-activity relation also for Ha and that binary com- 
ponents are stronger dependent on period than single stars. 
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6.5 An activity-Li abundance relation? 



In Fig. 16 1 and Fig. 16? the lithium abundance relative to 
hydrogen is shown as a function of effective temperature 
and as a function of rotational period, respectively. The sam- 
ple consist of our 18 singles and 42 binaries from Table |7] 
(the triple systems were omitted), 236 of the 408 CABS bi- 
naries that have a Li abundance, plus 75 out of the 109 single 



stars from paper I. Also included in Fig. 16 are the young 



solar analogs from Gaidos et al. (20001. Seven stars in the 



CABS catalog were in common with the more recent Li 
study by Maldonado et al. (2010 ) but only one of them had a 
reasonable lithium equivalent-width measure (HIP 78843), 
which we converted to a logarithmic abundance of 0.72 with 
Tefi=4790K and the NLTE tables in Pavlenko & Magazzu 



(1996 1. The same conversion was applied to the averaged 
equivalent-width measures for 78 single stars and 6 bina- 
ries from Lopez-Santiago et al. (2010) and 153 singles and 
24 binaries from Maldonado et al. (20 10) and added to the 
sample in Fig. 16 1. For these stars, only B — V colors are 
listed and Tcff 's are consequently obtained from the B — V 
calibration of Flower (1996) . Despite that B — V colors are 
affected by metallicity (see, e.g.. Gray |1994) l, and that at 
least Maldonado et al. ( |2010| l obtained metallicities from 
their spectra for most of their targets, the typical metallicity 
range for our stars is at most ±0.3 dex with the one or other 
outlier. This would affect the temperature scale at most by 
±50 K and therefore, we decided not to correct for metallic- 
ity simply because metallicities were not known for all stars 
in our samples. 

6.5.1 The lithium-temperature dependency 



Figure [T6p shows a generally increasing Li abundance for 
binaries with increasing effective temperature (and there- 
fore likely with increasing mass). Such a result was found 
earlier for single stars in various open clusters, most note- 
worthy in the Hyades, Pleiades, Praesepe, M34, and M67 
(see, e.g., Jones et al. |1997[ ). Our binary sample shows in- 
creasing dispersion of Li abundance with increasing tem- 
perature, opposite to what was found for Pleiades stars (e.g. 



Soderblom et al. 1993a i but in agreement with the Hyades 
FGK stars (e.g. Thorburn et al.'1993). It reaches up to three 
orders of magnitude for binaries with Tcff «5000 K, i.e. 
15-20(7 of the observational uncertainty, and still two orders 
of magnitude between 4500-5000 K. The dispersion is still 
one order of magnitude at the cool end at Tcff < 4500 K. 
For the binaries hotter than 5000 K, the maximum observed 
Li abundances are primordial (sa 3.5) and the dispersion 
decreases due to the generally higher abundances with tem- 
perature. Part of the explanation for such a drastic disper- 
sion is certainly that our binary sample is a mix of ages, 
most likely ranging from ZAMS binaries to maybe 10 Gyr 
systems. This is evidenced by the fact that our single-star 
sample in Fig. [T6^ shows a qualitatively comparable behav- 
ior but the dispersion never reaches much above two orders 
of magnitude. 



Binaries have larger systematic errors of the equivalent- 
width measures than singles, and that certainly also adds 
to the dispersion. It is particularly due to the double-lined 
nature for some binaries and their relative continuum uncer- 
tainties and associated larger errors for Tcff. However, we 
expect equivalent- width errors never larger than 10% for 
stars with equivalent widths > 50 mA and at most ss30% 
for equivalent widths < 20 mA while effective temperatures 
are good to within, say, ±50-70 K for singles and SB Is, and 
±100-150 K for SB2s. We also note that we use Li detec- 
tions and never just upper limits, and use temperatures de- 
termined from spectrum synthesis rather than photometry. 
Our Li abundances are therefore generally good to within 
±0.2 dex. Maybe other - extrinsic - factors may affect the 
Li abundance to an unknown amount but were deemed un- 
reasonable by Soderblom et al. (1993a); c/o their discussion. 
We conclude that a significant part of the Li-abundance dis- 
persion in active binaries must be intrinsic. 

6.5.2 The lithium-rotation dependency 

If the depletion mechanism(s) were dependent on stellar ro- 
tation, the expectation is that tidally locked binaries would 
exhibit lower surface abundances than single stars (of same 
mass and age). If independent of rotation, one could expect 
higher Li abundance in binaries because of less depletion. 



There is evidence in Fig. 16 i that binaries exhibit on aver- 
age less surface lithium than singles. Log-linear fits to the 



data of the form log n(Li) 
average relations; 



aTpfj + c provide the following 



log n(Li) =aT^ 



cff 



0.0010405 _ ^ 



for binaries, (13) 



logn(Li) = c '"-'-'^ - d for singles, (14) 

where a, b, c, d are the constants 58.29306 10^, 58.8000 10^, 
7.66085 10^, 8.1658010^ respectively. The difference be- 
tween the two fits is significant for our sample and amounts 
to wO.25 dex at 5500 K. However, the fits should not to be 
taken too literally and act just like a guideline. 

The general picture is that the Li depletion is related 
to the star's angular momentum loss, itself related to inter- 
nal flows and shear instabilities, mixing, as well as external 
winds (e.g. Pinsonneault et al. |1990| l. A crystal-clear picture 
still has to emerge but that magnetic fields and winds play 
a major role has been accepted widely and new numerical 
simulations are going in this direction (e.g. Matt et al. 2012[ l. 
The rotational angular momentum must be thought of being 
coupled to the orbital angular momentum in binaries where 
the depletion mechanism(s) can chew on the orbital energy. 
There is some observational evidence that tidally locked bi- 
nary systems in the Hyades and M67 have larger lithium 
abundances than single stars though (see discussion in Bar- 



rado y Navascues et al. 2001 1 but, so far, did not amount to 
a conclusive (quantitative) relation. 

In Fig. [T6] 3, we plot Li abundance versus rotational pe- 
riod. Rotation periods range between 0.5 d and 200 d. Be- 
cause the vast majority of the binaries in our sample cluster 
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Fig. 16 Average lithium abundances a versus effective temperature, b versus (logarithmic) rotation period. Filled symbols 
refer to binaries, open symbols to single stars. The symbols and the sample sources are identified in the inserts. Active 
binaries show generally lower Li abundances than comparable singles. Our sample of binaries and single (field) stars 
exhibits a dispersion of up to 2-3 orders of magnitude in the full range 4500-6000 K that peaks at «5000 K. When plotted 
versus (logarithmic) rotation period, the dispersion is even 3^ orders of magnitude. 



below or near periods of w 1 d, we had to choose a logarith- 
mic presentation for the period axis as well. Nevertheless, 
the log-log plot shows a clear trend of higher Li abundance 
with faster rotation, for binaries as well as for singles, main 
sequence or evolved. Again, as in the Li-Toff plane, the dis- 
persion is very large, peaking at 3^ orders of magnitude 
at around a rotation period of 10 d. This effectively camou- 
flages a clear Li-Prot relation, in particular because we do 
not know ages for our systems as one does for open clusters. 
As an averaged guideline, we fit linears to the binary data in 
the log-log plane and obtain 

logn(Li) =a logProt + b for binaries, (15) 

where a, h are the constants —0.60067 and 1.84783, respec- 
tively. No fit is given for the single stars because our sample 
is not complete due to a general lack of rotation periods. 
Jones et al. ( |1997[ ), together with parallel and subsequent 
work by several other groups, presented clear evidence that 
the (single) rapid rotators in the Pleiades (age «100 Myr) 
have generally also higher Li abundances than their slowly 
rotating cousins, at least for the temperature range 4400- 
5400 K. A similar trend appeared in the data for M34 (age 
w300 Gyrs) but only for stars of w4500 K, while the Hyades 
(625 Myr) shows no rapid rotators anymore and thus no Li 
relation with rotation. 

Several mechanisms were suggested to explain the high 
Li abundances in rapidly rotating, active stars. Among them 
the production of Li during coronal flares through spalla- 
tion reactions, or through an underestimated surface tem- 
perature due to large amounts of cool starspots (see, e.g., 
Pallavicini et al. 11993) or simply through the explanation 
that active stars are generally being young (Randich et al. 
|1993[ ). Another interesting explanation for tidally coupled 



binary components is through a lack of rotationally induced 
mixing due to the little magnetic braking compared to sin- 



gle stars (Pinsonneault et al. 1990 1. Bohm-Vitense (2004 1 
noted that for giants with Tcff less than sa6300 K the re- 
duced lithium surface abundance appears related with the 
decreased surface rotation. She attributes this to effects from 
deep mixing. Just recently, Takeda et al. (2010) announced 
evidence for a (positive) correlation of Li abundance with 
rotational velocity in solar-analog stars. A depletion of two 
orders of magnitude in Li abundance is seen in the range of 
V sin i from 4 km s^^ to 1 km s^^ according to their Fig. 5. 
It cafls for a very efficient envelope mixing process, possi- 
bly also related to the exterior angular momentum carried 
by planets. 

Accretion of already a few Earth masses during early 
main-sequence phase would bring up the ^Li abundance 
in these stars and a proper amount of mixing may prevent 
^Li from destruction (Murray et al. 2001 1. Such planet en- 
gulfing would be in agreement with the observation of in- 
creased Li abundance in rapidly-rotating active stars but re- 
quires the detection of ^Li in active stars. Such a detec- 
tion was reported by Israelian et al. (2001 ) for the planet- 
hosting solar-type star HD 82943 or the active subgiant bi- 



nary HD 12353 1 by Strassmeier et al. ( 201 1 1 (but see Cayrel 



et al. [2007 1 for a warning of such detections). 

7 Summary and conclusions 

We present and analyzed high-resolution spectra from the 
STELLA robotic telescope and determined orbital elements 
for 45 binaries in the period range 0.29-6450 d, many for 
the first time but all with substantially increased precision 
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and likely also increased accuracy. Additional photometric 
monitoring of many of these stars with our APTs allowed 
the determination of precise rotation periods. 

We summarize our findings as follows. 

- STELLA/SES achieves a radial-velocity precision for 
an individual radial- velocity observation of »30 ms~^ 
without an iodine cell or a simultaneously recorded 
comparison spectrum. It enables rms values of as small 
as «40 m s~^ for fits of orbital elements. For one of the 
favorable cases in this paper (the SB2 binary HIP 77210 
with Porb=9.9 d) the minimum primary mass is precise 
to 0.12% and the secondary mass to even 0.07% . 

- The radial-velocity zero-point difference be- 
tween STELLA and CORAVEL measurements is 
H-0.503kms-^ 

- The bulk (74%) of rapidly-rotating active stars in bina- 
ries are synchronized and in circular orbits. 

- About 26% (61 targets) of the binaries in our full sam- 
ple are asynchronous rotators. About half of them have 
^rot > ^orb and of these, all but two have e > 0. 
As rotational synchronization is predicted to occur be- 
fore orbital circularization active binaries must have 
gone through an extra spin-down phase besides tidal 
dissipation and we suggest this to be due to a mag- 
netically channeled wind with its subsequent braking 
torque. Such braking depends on the stellar magnetic- 
field geometry and would result in very different brak- 
ing efficiencies for stars of otherwise nearly identical 
physical parameters. 

- We find a steep increase of Prot for lower T^s in both 
single stars and binary components. This verifies ear- 
lier claims for active single stars with known cycles and 
shows that it is also the case in our active binaries. A 
functional dependence of Prot oc is suggested, 
where X = 7.2 for single stars and X = 7.0 for bi- 
nary stars. The difference in X is not significant though. 

- A relation between Ha-core flux and rotational period 
is evident for both single stars and binaries. Power- 
law fits suggest Rua oc Prot binary stars and 
-Rho oc Pfot '^^ for single stars, where R is Ha-core flux 
normalized to (tT^J'jj . The difference of the exponents is 
only weakly significant on the 3a level. 

- Our data suggest that Li abundances of active binary 
components also increase with effective temperatures, 
as known for single stars. However, binaries appear to 
have on average 0.25 dex lower Li abundances than sin- 
gle stars of same effective temperature, which is what 
we expect if the depletion mechanism is rotationally en- 
hanced in binary stars. The dispersion for binaries ap- 
pears to peak around Tea ~5000 K, amounting to three 
orders of magnitude. This dispersion decreases towards 
both ends of the temperature range between 4500 K to 
6500 K. Our single-star sample follows this trend as 
well but with an order of magnitude smaller dispersion. 
We can not separate age effects in our sample but at- 
tribute a significant fraction of the dispersion to it. 



- Rotational dependency of Li abundance in binaries is 
also evident and suggests a law of form logn(Li) oc 
—0.6 log Prot as a guideline. However, a dispersion of 
up to 3-4 orders of magnitude explains why no quantita- 
tive relations were found so far. The dispersion appears 
largest for binaries with periods shorter than «10 d and 
still amounts to 2-3 orders of magnitude for periods 
larger than «10 d. 
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